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Abstract. We analyze a class of weakly diflferentiable vector fields F: — > with the prop- 
erty that F G L°° and divF is a Radon measure. These fields are called bounded divergence- 
measure fields. The primary focus of our investigation is to introduce a suitable notion of the 
normal trace of any divergence-measure field F over the boundary of an arbitrary set of finite 
perimeter, which ensures the validity of the Gauss-Green theorem. To achieve this, we first de- 
velop an alternative way to establish the Gauss-Green theorem for any smooth bounded set with 
F S -L°°. Then we establish a fundamental approximation theorem which states that, given a 
Radon measure fj, that is absolutely continuous with respect to 7-{^~i on M.^ , any set of finite 
perimeter can be approximated by a family of sets with smooth boundary essentially from the 
measure- theoretic interior of the set with respect to the measure \\fJ,\\. We employ this approxi- 
mation theorem to derive the normal trace of F on the boundary of any set of finite perimeter, 
E, as the limit of the normal traces of F on the boundaries of the approximate sets with smooth 
boundary, so that the Gauss-Green theorem for F holds on E. With these results, we analyze the 
Cauchy fluxes that are bounded by a Radon measure over any oriented surface (i.e. an (A^ — 1)- 
dimensional surface that is a part of the boundary of a set of finite perimeter) and thereby develop 
a general mathematical formulation of the physical principle of balance law through the Cauchy 
flux. Finally, we apply this framework to the derivation of systems of balance laws with measure- 
valued source terms from the formulation of balance law. This framework also allows the recovery 
of Cauchy entropy fluxes through the Lax entropy inequality for entropy solutions of hyperbolic 
conservation laws. 



1. Introduction 

In this paper we analyze a class of weakly differentiable vector fields F : M.^ with the 

property that F G L°° and divF is a Radon measure fi with finite total variation (i.e. a totally 
finite signed measure). These fields are called bounded divergence-measure fields, and the class is 
denoted by VA4°°. The primary focus of our investigation is to introduce a suitable notion of 
the normal trace of any divergence-measure field over the boundary of an arbitrary set of finite 
perimeter to obtain a general version of the Gauss-Green theorem. Clearly, this investigation is 
closely related to the theory of BV functions in R^; in fact, it would be completely subsumed by 
the BV theory if the fields were of the form F = (i^i, F2, . . . , Fm) with each Fk G BV{R^), since 
then divF = J2k=i ff^ ^ Radon measure fi with finite total variation (cf. [Zic89]). However, in 
general, the condition divF ~ allows for cancellation, which thus makes the problem more difficult 
and accordingly more important for applications (see §9-§ll). For the Gauss-Green theorem in the 
BV setting, we refer to Burago-Maz'ja [BM69], Volpert [Vol67a], and the references therein. The 
Gauss-Green theorem for Lipschitz vector fields over sets of finite perimeter was first obtained by 
DeGiorgi [DG61a, DG61b] and Federer [Fed45, Fed58]. Also see Evans-Gariepy [EG92], Lin- Wang 
[LY02], and Simon [Sim83]. 
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Some earlier efforts were made on generalizing the Gauss-Green theorem for some special situ- 
ations of divergence- measure fields, and relevant results can be found in Anzellotti [Anz83] for an 
abstract formulation when F E L°° over a set with boundary and Zicmer [Zie83] for a related 
problem for div F € L°°; also see [AFPOO, BC84, BF91, JN90, JN94b, JNm95, Non95, Pfe05a, 
DPP04a, PfeOl, LRX07]. In Chen-Frid [CF99a, CF03], an explicit way to formulate the suitable 
normal traces over Lipschitz dcformablc surfaces was first observed for F G VA4°°. In particular, it 
was proved in [CF99a, CF03] that the normal trace over a Lipschitz deformable surface, oriented by 
the unit normal vector v, is determined completely by the neighborhood information from the posi- 
tive side of the surface pointed by v and is independent of the information from the other side. This 
is the primary motivation for our further investigation on divergence-measure fields. Chen- Torres 
[CT05] were the first to obtain the normal trace for any bounded divergence-measure field over a set 
of finite perimeter. E, and the corresponding Gauss-Green theorem. One of the main results in this 
paper is to obtain this normal trace as the limit of the normal traces over the smooth boundaries 
that approximate the reduced boundary d*E of E. In particular, the normal trace is determined 
completely by the neighborhood information essentially from the measure-theoretic interior of the 
set (see Theorem 21.1), so that the Gauss-Green theorem holds for any set of finite perimeter. 

We recall a very general approach, initiated by Fuglede [Fug60] , in which the following result was 
established: If F e Lp(IR-'^; IR-'^), 1 < p < oo, is a vector field with divF = fi, where fj. is a signed 
Radon measure, then 

(2.1) / divF := m(£) = - / F{y) • u{y) dH^-'iy) 

for "almost all" sets of finite perimeter, E^ where Ti^^^ is the {N — l)-dimensional HausdorfF 
measure. The term "almost all" is expressed in terms of a condition that resembles "extremal 
length", a concept used in complex analysis and potential theory (cf. [Zie67, Zie69, Zie70, HcOT]). 
One way of summarizing our work in this paper is to say that we wish to extend Fugiede's result 
so that (2.1) holds for every set E of finite perimeter. Of course, this requires a suitable notion of 
the normal trace of F on d* E. This is really the crux of the problem as F, being only measurable, 
cannot be re-defined on an arbitrary set of dimension N — 1. To achieve our goal, we first establish 
a fundamental approximation theorem which states that, given a Radon measure n on such 
that n << TC'^~^, any set of finite perimeter can be approximated by a family of sets with smooth 
boundary essentially from the measure-theoretic interior of the set with respect to the measure 
(e.g. fj, = divF). Then we employ this approximation theorem to derive the normal trace of F on 
the boundary of any set of finite perimeter as the limit of the normal traces of F on the smooth 
boundaries of the approximate sets and establish the Gauss-Green theorem for F which holds for 
an arbitrary set of finite perimeter. 

With these results on divergence-measure fields and sets of finite perimeter, we analyze the Cauchy 
flux that is bounded by a nonnegative Radon measure a over an oriented surface (i.e. an {N — 1)- 
dimensional surface that is a part of the boundary of a set of finite perimeter) and develop a general 
mathematical formulation of the physical principle of balance law through the Cauchy flux. In the 
classical setting of the physical principle of balance law, Cauchy [Cal823, Cal827] first discovered 
that the flux density is necessarily a linear function of the interior normal (equivalcntly, the exterior 
normal) under the assumption that the flux density through a surface depends on the surface solely 
through the normal at that point. It was shown in Noll [Nol59] that Cauchy's assumption follows 
from the balance law. Ziemer [Zie83] provided a first formulation of the balance law for the flux 
function F £ L°° with divF € L°° at the level of generality with sets of flnite perimeter. Also see 
Dafermos [Daf98] and Gurtin-Martins [GW67, GW71]. One of the new features in our formulation 
is to allow the presence of exceptional surfaces, "shock waves", across which the Cauchy flux has 
a jump. When the Radon measure a reduces to the Af-dimensional Lebesgue measure £^ , the 
formulation reduces to Ziemer 's formulation in [Zie83], which shows its consistency with the classical 
setting. We flrst show that, for a Cauchy flux bounded by a measure a, there exists a bounded 
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divergence- measure field F : R — > R , defined £ -a.c., sucli tliat 



Js 



for almost any oriented surface S, oriented by the interior unit normal ly. Then we employ our 
results on divergence-measure fields to recover the values of the Cauchy flux on the exceptional 
surfaces directly via the vector field F. The value as the normal trace of F on the exceptional 
surface is the unique limit of the normal traces of F on the non-exceptional surfaces essentially from 
the positive side of the exceptional surface pointed by u. Finally, we apply this general framework to 
the derivation of systems of balance laws with measure- valued source terms from the mathematical 
formulation of balance law. Wc also apply the framework to the recovery of Cauchy entropy fluxes 
through the Lax entropy inequality for entropy solutions of hyperbolic conservation laws by capturing 
entropy dissipation. 

We observe the recent important work by Bourgain-Brezis [BB02] and De Pauw-Pfeffer [DPP04a] 
for the following problem with different view of point: Find a continuous vector field to the 

divergence-measure equation: 

(3.1) divF = fi inn, 

for a given Radon measure ii. In the case dji = f dx where / G the existence of a solution 

F to (3.1) follows form the closed-range or Hahn-Banach theorem as shown in [BB02]. It is proved 
in [DPP04a] that equation (3.1) has a continuous weak solution if and only if is a strong charge, 
i.e., given e > and a compact set if C fi, there is 6* > such that 



for any smooth function (j) compactly supported on K. 

The organization of this paper is as follows. In §2, we first recall some properties of Radon 
measures, sets of finite perimeter, and related BV functions, and then we introduce the notion of an 
oriented surface and develop some basic properties of divergence-measure fields. In §3, we develop 
an alternative way to obtain the Gauss-Green formula for a bounded divergence-measure field over 
any smooth boundary by a technique, which motivates our further development for the general case. 
In §4, we establish a fundamental approximation theorem which states that, given a Radon measure 
H on such that fi « H^~^, any set of finite perimeter can be approximated by a sequence 
of sets with smooth boundary essentially from the interior of the set with respect to the measure 
In §5, we introduce the normal trace of a divergence-measure field F on the boundary dE 
of any set of finite perimeter as the limit of the normal traces of F on the smooth surfaces that 
approximate dE essentially from the measure-theoretic interior of E with respect to the measure 
||diVF||, constructed in §4, and then we establish the corresponding Gauss-Green theorem. In §6-§7, 
we further analyze properties of divergence-measure fields, especially showing the representation of 
the divergence- measures of jump sets via the normal traces and the consistency of our normal traces 
with the classical traces (i.e. values) when the vector field is continuous. In §8, we first show that, if 
the set of finite perimeter, E, satisfies (29.1) (which is similar to Lewis's "uniformly fat" condition 
in potential theory [Lew88a]), there exists a one-sided approximation to E, and we then show that 
an open set of finite perimeter is an extension domain for any bounded divergence-measure field. In 
§9, we first introduce a class of Cauchy fluxes that allow the presence of these exceptional surfaces or 
"shock waves" , and we then prove that such a Cauchy fltix induces a bounded divergence-measure 
(vector) field F so that the Cauchy flux over every oriented surface with finite perimeter can be 
recovered through F via the normal trace over the oriented surface. 

In §10, we apply the results established in §3-§9 to the mathematical formulation of the physical 
principle of balance law and the rigorous derivation of systems of balance laws with measure- valued 
source terms from that formulation. Finally, in §11, we apply our results to the recovery of Cauchy 
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entropy fliix through the Lax entropy inequahty for entropy solutions of hyperbohc conservation 
laws by capturing entropy dissipation. 

2. Radon Measures, Sets of finite perimeter, and Divergence-Measure Fields 

In this section we first recall some properties of Radon measures, sets of finite perimeter, and 
related BV functions (also cf. [AFPOO, EG92, Fed69, Giu77, Zie89]). We then introduce the notion 
of oriented surfaces and develop some basic properties of divergence-measure fields. For the sake 
of completeness, we start with some basic notions and definitions. First, denote by H^^ the M- 
dimensional Hausdorff measure in for M < N, and by the Lebesgue measure in (recall 
that £^ = H^). For any -measurable set E C , we denote \E\ as the £''^-Lebesgue measure 
of the set E and dE as its topological boundary. Also, we denote B{x, r) as the closed ball of radius 
r and center at x. The symmetric difference of sets is denoted by 

AAB := (A \ B) U (B \ A). 

Let Q, C be open. We denote by S g that the closure of E is compact and contained in O. Let 
Cc(ri) be the space of compactly supported continuous functions on Q. with ||</5||o.q := sup{|(^(j/)| : 
2/ e fi}. A Radon measure in O is a regular Borel measure whose total variation on each compact 
set K Vl \s finite, i.e. ||/i||(if) < oo. The space of Radon measures supported on an open set Q. 
is denoted by M.{Vl). Any Radon measure /i can be decomposed into the difference of two positive 
Radon measures /i = /i^ — /^^; the total variation of /U is = /x"*" -|- Equivalently, if /U is a 
signed Radon measure on O, the total variation of fi on any bounded open set S C f2 is equal to 

Uoo 
^d^i:^eC,{B), < l} =sup{^Ki?,)|} 

i=0 

oo 

where the second supremum is taken over all pairwise disjoint Borel sets Bi with B = {J Bi. Since 

i=l 

the space of Radon measures can be identified with the dual of Cc(^^)i we may consider a Radon 
measure as a linear functional on Cc{fl), written as 

(4.2) := I (pdfi for each 9? e Cc(n). 

Jn 

We recall the familiar weak*-topology on M{Sl) which, when restricted to a sequence {nk}, yields 

IJ,k IJ- in M{0,), 
that is, fik converges to fi in the weak* topology if and only if 

(4.3) Mfe(<^) ~^ /^(v') for each <^ e Cc(ri). 

The space L''{ihfi,), I < p < oo, denotes all the functions / with the property that |/|^ is /i- 
intcgrablc. The conjugate of j; is q := p/{p — 1). The norm of / on a set E with integration 
taken with respect to a measure n is denoted by ^. In the event /x is Lebesgue measure, we 

will simply write H/llp.j^. 

4.1. Theorem (Uniform boundedness principle). Let X be a Banach space. IfTk is a sequence of 
linear functionals on X which converge weak* to T. Then 

limsup ||Tfe|| < 00. 

fc— *oo 

This theorem implies the following corresponding result for Radon measures. 

4.2. Corollary. Let /ifc he a sequence of Radon measures that converge to fi in the weak* topology. 
Then 

limsup WnkW < oo. 
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Next, wc quote a familiar result concerning weak*-convergence. 

4.3. Lemma. Let jj, be a Radon measure on fl and let Hk be a sequence of Radon measures converging 
weak* to fi. Then 

(i) If A c ft is any open set and fik o-re positive Radon measures, 

li{A) < liminf Hk{A); 

fc— *oo 

(ii) If K cfl is any compact set and Hk are positive Radon measures, 

IJ-i^) > limsup/ife(ii'); 

fe— *oo 

(iii) // ||/Ufe|| a, then < a. In addition, if E satisfies a{dE) = 0, then 

tx{E) = lim f,k{E). 

More generally, if f is a bounded Borel function with compact support in Vt such that the set of its 
discontinuity points is a -negligible, then 

lim / fdnk= / / d^. 



5.1. Definition. For every a e [0, 1] and every -measurable set E c M^, define 
(5.1) £;«:={yeM^ : D(^,y)=a}, 

where 

Then E" is the set of all points with density a. We define the measure-theoretic boundary of E, 
d'^E, as 

(5.3) d"'E:=R^\{E°UE^). 

5.2. Definition. A function /: — > R is called a function of bounded variation if each partial 
derivative of f is a totally finite signed, Radon measure. Notationally, we write f G BV{^V). Let 
E Vl be an -measurable subset. We say thai E is a set of finite perimeter if Ti.^ ~^ {d"^ E) < oo. 
Equivalently, E is of finite perimeter if G BV{^V). Consequently, if E is a set of finite perimeter, 
then VXg is a (vector-valued) Radon measure whose total variation is denoted by HVX^jj. 

5.3. Definition. Let E <^ Cl be a set of finite perimeter. The reduced boundary of E, denoted as 
d*E, is the set of all points y € fl such that 

(i) \\VxJiBiy,r))>Oforallr>0; 

(ii) The limit u^{y) := lim^^o ||vX^{|tBfa,r)) 
Then, for Ti^-^-a.e. y G d*E, 

r^o a{N-l)r^-^ 

where a{N — 1) is the Lebesgue measure of the unit ball in and the generalized gradient ofX^ 

enjoys the following basic relationship with H^~^: 

(5.4) \\VxJ=n^-^Ld*E. 

The unit vector, (y), is called the measure-theoretic interior unit normal to E aty (we sometimes 
write V instead of ve for notational simplicity). Also, we recall that the reduced boundary, d*E, is 
an {N — l)-rectifiable set which implies that there exists a countable family of -manifolds of 
dimension N — 1 and a set Af ofH^~^ measure zero such that 

oo 

(5.5) d*Ec{[j Mk)[jJ^- 

fe=i 
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In view of the following, we see that u = is, aptly named because u is the interior unit normal 
to E provided that E (in the limit and in measure) lies in the appropriate half-space determined by 
the hyperplane orthogonal to v; that is, u is the interior unit normal to at a; provided that 

D{{y : {y - x) ■ u > 0,y ^ E}U {y: {y - x) ■ u < 0,y e E}, y) = 0. 

We will refer to the sets E'^ and E^ as the measure-theoretic exterior and interior of E. Wc note 
that, in general, the sets and E^ do not coincide with the topological exterior and interior of 
the set E. The sets E'^ and E^ also motivate the definition of measure-theoretic boundary. Indeed, 
for any set E C R^, the definitions imply that ^ E^ U d*E U E° U Af where n'^-^{J\f) = 0. If 
we define a set E to be "open" if E is both measurable and D{E,x) = 1 for all x £ E, then this 
concept of openness defines a topology, called the density topology. It is an interesting exercise to 
prove that the open sets are closed under arbitrary unions; the crux of the problem is to prove that 
the arbitrary union is, in fact, measurable. This topology is significant because it is the smallest 
topology (the one with the smallest number of open sets) for which the approximately continuous 
functions are continuous [GNN61]. 

6.1. Remark. If E is a set of finite perimeter, then clearly 

(6.1) d*E cEi cd^'E, n^-'^{d'^E\d*E)=0. 

The following result, which is easily verified (although tedious), will be needed in the sequel. 

6.2. Lemma. If A,B (^Q. are sets of finite perimeter, then 

a'"(A r\B) = (a™ A n s) u (A n a™B) u (a™ A n a^s) . 

6.3. Definition. Let p S C^{M.^) be a standard symmetric mollifying kernel; that is, p is a non- 
negative function with support in the unit ball and satisfies \\p\\i.^N = 1. With u S L^{R.'^), we set 
Ue := u* Ps, where the sequence Ps{y) '■— prp(f ) forms a mollifier. 

Recall the following properties of mollification (cf. [Zie89]): 

6.4. Lemma. (i) If u & L^„^(K^), then, for every £ > 0, e C~(M^) and D'^{p^ * u) = 

{D°'ps) * u for each multi-index a; 
(ii) With k := 1/ eu and Sk ^ 0, Wfe(.x) u{x) whenever x is a Lebesgue point of u. In particular, 
if u is continuous, then converges uniformly to u on compact subsets of M.^ . 

When u is taken as for a set of finite perimeter, E, this result can be considerably strengthened. 

6.5. Lemma. Ifuk is the mollification ofx^ for a set of finite perimeter, E, then the following hold: 

(i) Uk e C°°(R^); 

(ii) There is a set M with 'H^~^{M) = and a function e BV such that, for all y ^ Af, 
Uk{y) — *■ Uj^{y) as fc ^ 00 and 



1 y&E\ 
\ y&d*E, 
yGE°; 



(iii) Vuk Vu^ inM{m.^); 

(iv) ||Vufe|| {U) — > ||Vue|| (f/) ask ^ oo, for any open set U with ||Vue|| {dU) = 0; 

(v) VX^ = Vw^. 

Proof. Only (iii) requires a proof, since (i), (ii), and (iv) are the results from the standard BV 
theory and (v) is immediate from the definitions and the fact that Uj^ = almost everywhere. As 
for (iii), since Uk in L^(M^), then Uk when considered as distributions, which implies 

that Vufc as distributions and consequently as measures since Vufe, Vu^ G M.{R^). □ 
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6.6. Remark. More generally, functions in the spaces BV{R^) and W^'P{R^), 1 < p < oo, have 
precise representatives; that is, if u € BVi^^), then there is a function u* e BV{U-^) such that u 
and u* are equal a.e. and that the mollification sequence of u, Uk, converges to u* at all points except 
those that belong to an exceptional set E with TL^~^{E) = 0. However, this is not the same as saying 
that u has a Lebesgue point, which is slightly stronger. A similar statement is true for functions in 
the Sobolev space W^'^iM.^), 1 < p < oo, except that the exceptional set E has jp-capacity zero, see 
Definition 8.2 below. As we will see, the ji-capacity vanishes precisely on sets ofTi,^~^ measure 
zero. Thus, we can say that functions in the spaces BV and W^'^ have precise representatives that 
are defined, respectively, 71 and 7p almost everywhere. 

The next result affirms the notion that the mollification is generally a norm reducing operation. 

7.1. Lemma. Let E be a set of finite perimeter and let denote the mollification of X^^. Then 

\Wuk\\^ < \\VXJ\. 

Proof. For any / G BV{R^), consider the convolutions fe{y) = /rn Peiv — x)f{x) dx. Using V/e = 
* V/ and fe e C°°(]R^), we obtain 

^fe{y)= / pe{y - x) dm{x) 

where m := V/ is the measure. Thus, we have 

|V/e(y)|< / p,{y-x)d\\m\\{x). 
In particular, when / = and f^^ = Uk with Sk = 1/fc, then m = Vxe and 
\^uk{y)\< I pe^{y-x)d\\m\\{x) for all y G M^. 

That is, 

/ \^Uk{y)\dy< / / pe{y - x) d\\m\\{x) dy 

= / / pe{y-x)dyd\\m\\{x)<\\mm''). □ 

We recall that the BV space, the space of functions of bounded variation, in fact represents 

equivalence classes of functions so that, when a function in a class is changed on a set of C^- 
measure zero, it remains in this class. The same is true for sets of finite perimeter because, by 
definition, the characteristic function of a set of finite perimeter, E, is a function of bounded 
variation. Thus, it follows that E may be altered by a set of /3^-measure zero and still determine 
the same essential boundary d"^E. Throughout, we will choose a preferred representative for E and 
thereby adopt the following convention. 

7.2. Definition. E := {y : D{E,y) = 1} U d'^E. 

7.3. Definition. A vector field F G 1 < p < 00, is called a divergence-measure field, 
written as F £ 'DMp{^), if p := divF is a (signed) Radon measure with finite total variation on 
in the sense of distributions. Thus, for (p G C^{n), we have 

p{(p) := divF{(f) = - / F -Vipdy. 

Jn 

The total variation of p is a positive measure which, for any open set W , is defined as 
\\p\\ {W) : = supipicp) : M,.^^ <l,^e C^{W)} 

= sup {J^F-Vipdy : M^.^ < 1, ^ e (W)} . 

A vector field F G Lf^^i^l) means that, for any K Q., F G VMp{K). 
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8.0. Definition. Let F G T>Mp{^1), 1 < p < oo. For an arbitrary measurable set i? g fJ, the trace 
of the normal component of F on dE is a functional defined by 



(8.1) 



{TF),,{^) = [ Vip-Fdy+ [ ifdii 
Je Je 



for all test functions ip £ C^{Q). Clearly, (TF)gj^ is a distribution defined on fl. Note that this 
definition assumes only that the set E is measurable. Later, we will provide an alternative definition 
when E is a set of finite perimeter (see Theorem 25.1). 

8.1. Proposition. Let E Q be an open set. Then spt{{TF)gj^) c dE. That is, if tjj and (p are 

test functions in D(f2) with ij) = on dE, then (TF)gg(^) = {TF)gj^{ip). 

Proof. If the support were not contained in dE, there would be a point xq ^ dE with xq G 
spt{{TF)g^) n E. This imphes that, for each open set U containing Xq, there exists a test function 
ip G C^{U nE) such that {TF)g^{ip) ^ 0. Choose U so that [/ c \ dE. Let F^ denote the 
molUfication of F (see Lemma 6.4). Then, since spt{Fs(p) <e E, 



0= / diviF,ip)dy = {TF,),,{ip) 

J E 

= Fe-V<p+ ipdivF^dy 
Je Je 

^ / F-Vipdy+ [ ipdn={TF)g^{ip)i^Q, 
J e Je 



' E 

where we used dE n spt((p) = in the Umit. Thus, we arrive at our desired contradiction. □ 

8.2. Definition. For I < p < N , the p-capacity of an arbitrary set A g is defined as 
(8.2) 7p(A):=inf{^|V<prrfy}, 

where the infimum is taken over all test functions ip G (Q) that are identically one in a neighbor- 
hood of A. It is well known (cf [FZ73a]j that jp{A) = for 1 < p < N implies that H^-p+^{A) = 
for each £ > and that, conversely, if TI^~p{A) < oo, then jp{A) = 0. In view of Remark 6.6 and 
Lemma 9.2, it is easy to verify that the class of competing functions in (8.2) can be enlarged to the 
Sobolev space W'^'P{Q). 

8.3. Remark. The case of p = 1 requires special consideration. In 1957, Fleming conjectured that 
7i(A) = if and only ifTi,^~^{A) = 0. This was settled in the affirmative by Gustin [Gus60] who 
proved the boxing inequality, from which Fleming's conjecture easily follows (cf. [Fle60]^. 

The next result is basic (cf. [Zie89, Zie69, ZieTO]). 

8.4. Proposition. Let 7p be the p-capacity defined as in the previous definition. Then 
(i) // Ek C is a sequence of arbitrary sets, then 

7p (lim inf Ek ) < lim inf jpiEk); 

k — *oo k — ^oo 



(ii) If El C E2 C ■ ■ ■ are arbitrary sets, then 

7p ( U -^fe ) = lim IpiEk); 

(iii) If Ki D K2 D ■ ■ ■ are compact sets, then 

{ f] ^k] = lim lp{Kk); 
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(iv) If {Ek} is a sequence of Borel sets, then 

7p I U -Sfe ) < lim JpiEk); 
Vfe=l / fc^oo 

(v) If A c M.^ is a Suslin set, then 

sup{7p(^r) : CA} = inf{7p(f/) : f/°P«" D A}. 

Any set function, 7, satisfying conditions (i)-(iv) is called a true capacity in the sense of Choquet 
and a set A satisfying condition (v) is said to he 7-capacitable. 

9.1. Remark. One of the main reasons for studying the capacity is its important role in the develop- 
ment of Sobolev theory. It was first shown in [FZ73a] that every function u G W^'P{fl) has a Lehesgue 
point jp-a.e. In particular, in view of Remark 8.3, this implies that a function u e W^^'^(f2) has a 

Lehesgue point everywhere except for an exceptional set E with 'H^~^{E) = 0. In case u € BViVL), 
we have a slightly weaker statement than the corresponding one for u G W^'^{^): 

lim + u{y) dy = u{x) for 'H^~^-a.e. x Gil,. 

"■^O iB(x,r) 

It turns out that the Sobolev space is the perfect functional completion of the space C^(f2) relative 
to the p- capacity. See [AS56a] where the concept of perfect functional completion was initiated and 
developed. 

9.2. Lemma. If F e VMf„^{R^), 1 < p < 00, then \\divF\\ (B) = whenever B ^ VLis a Borel 
set with 'yq{B) = 0, q -.= p/ [p— 1). In particular, when p = 00 (i.e. q= 1), then ||djVF|| (B) = 0. 

Proof. Because of the inner regularity of 7, and condition (v) of Proposition 8.4, it suffices to show 
that /i(-ft') = for any compact set K C B, where /i := divF. Since Jq{K) = 0, then there exists a 
sequence of test functions (pk G C^{Q) such that 

(i) (/5fe = 1 on K; 

(ii) ||V^,||,^0; 

(iii) (pk{y) for ally Gfl except those in some set A C il with 7g(vl) = 0. 
To see that such a sequence (pk exists, we proceed as follows: 

Case 1: q > 1. Since ^q{K) = 0, we may choose Ui G C^(0) with m = 1 near K and 
||Vwi||g ~^ 0- Then Sobolev's inequality implies 

\\ui\l<C\\Vur\\^ iorr^Nq/N-q. 

Thus, it follows that, for a subsequence, Ui{x) for a.e. x. In fact, an application of Mazur's 
theorem shows that the space of convex combinations of W^'P{VL is strongly closed and therefore, 
for a suitable subsequence, the sequence of finite convex combinations of the Ui, say Vk := J^ajUj 
with Yli=i fli = 1, it follows that Vk{x) for all x except those in an exceptional set A where 
7g(A) = (see [FZ73a], p. 156). Taking (pi := Vi for this subsequence establishes a sequence 
satisfying conditions (i)-(iii) above. 

Case 2: q = 1 {p = 00). Then the argument is modified by considering the following variational 
problem: 

(9.1) rk,qiK):=M\^J IVukl'dy : UkGW''^{R^),Uk = lonlh, spt{uk) CWk}. 

Let Wk D Uk D Uk - ■ ■ D K with Uk open and 

(9.2) 7p(Wfe) = 7p(Wfc) <ek^O. 
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We will also assume to be nested: Wi D W2 • ■ • ■ For each k, let Wi~ be a smoothly bounded set 
to find that f] Wk = K. Consequently, we may take (fk to be a suitably small mollification of X„,^ 
to sec that conditions (i) (iii) above arc satisfied. 

To prove our lemma, it suffices to show that, if K is any compact set satisfying jq{K) = 0, then 
IJ,{K) = 0. Since 

(a) </3fc = 1 on K; 

(b) (fk^O 7q-a-e.; 

(c) llV^Pfell^.aw ^0; 

(d) IImII {W) = sup{/^ \F -Wipldy.ipe C^{W), \ip\ < 1} for any open set W <^n, 
it follows that 



liin ( ipkdiJ,+ ipk dji) = fi{K) + li: 

Jk Juk\K f"- 

< n{K) + lim / 



Uk\K 



d\ 



Uk\K 



< lim 


/ ipkdivF 


< lim / \F 




Jn 


fe^oo Jn 



< llFll^^lim ||V¥'fe||^< ||F||^.0, 



thus obtaining ii{K) — 0, as desired. 



□ 



10.1. Corollary. If F e VM^^jn) for 1< p < N with divF =: /j, and ifn^-'i{A) < 00 for A ^ ^ 

with q = p/{p — 1), then 7q(^) = and hence ||/i|| (A) = 0. 

10.2. Remark. If F € T>M''{R^), 1 < p < 00, with divF = p. Then, in view of the fact that (p is 
defined 'H^~''-a.e. and therefore fx-a.e., with ip G W^''^{U.^), it follows that the integral 

I ipdji 

is defined and is meaningful. 

10.3. Example (Chen-Frid [CF99a]). Denote U the open unit square in M? that has one of its sides 
contained in the line segment 

L:={y = (^1,2/2) : yi = 2/2} nac/. 
Define a field F -.m^XL^m? hy 

F{y) = F{yi,y2) 



sm 



) , - sm ( 



yi - 2/2 



2/1 - 2/2 



))• 



Clearly, F 6 -L°°(R^), and a simple calculation reveals that divF = in \_L. Then F belongs 
to VM.°° {R"^); but the field is singular on one side, L, of dU and therefore, F is undefined on dU ; it 
has no trace on dU in the classical sense. Note also that the points of L are all essential singularities 
of F because the following limit docs not exist: 

lim F{y) {ov y \ L, x € L, 
and therefore the normal trace of F on dU is given by 

lim / F{y) ■ u{y) dH^iy) = hm / divFdy = lim = 0, 

t *0 J ^ t s-O J j-j^ t *0 

where Ut := {y & U : d{y, dU) > t)}. Thus, we have shown the following: 

(i) F is an element of X>jV1°°(K^), while each component function of F is not in BV{R'^); 

(ii) F has an essential singularity at each point of L and therefore cannot be defined on L; 
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(iii) As we will see later, Theorem 25.1, F has a weak normal trace on L which is sufficient for the 
Gauss-Green theorem to hold. 

For more properties of the spaces VM^ of divergence-measure vector fields, see Chen-Prid [CF99a, 
CF03]. 

The following theorem provides a product rule for the case p = oo. For the sake of completeness, 
we will include its proof, which is slightly different from that given in [CF99a]. We denote by {g^} 
a sequence of C^{M.^) mollifications with the property that gk ^ g in L^{M.^) and such that 
llVgfcll ^ IIV.9II (cf. [Zie89], p.500). 

11.1. Theorem (Chen-Frid [CF99a]). Let F e VM°°(R^) and g G BV{R^). Then 



(11.1) div{gF) = g* divF + F -Vg, 



where F ■ Vg denotes the weak* limit of the measures F-Vgu and g* denotes the limit of the mollifiers 
ofg (cf [Zie89];. 

Proof. Let F^ be the mollification of F and set fi := divF. Since are smooth, the classical 
product rule yields 

(11-2) div{gkF,) = gk divF, + F, ■ Vgk- 

First, we note that divFg = {divF)^ = fJ-s M in jVI(M^) as £ ^ 0. Since gkF^ — > gkF in 
Lj^^{M.^) as e ^ 0, we obtain from (11.2) that, in the sense of distributions, 

(11-3) F-Vgk = div{gkF)-n{gk). 

Owing to the fact that F G L°°, we see that F ■ Vgu is a bounded sequence in L^(R^) and hence 
there is a subsequence such that F ■ Vgk converges weak * to some measure, denoted by F ■ Vg. 
Letting fc ^ 00 in (11.3) yields 



(11.4) F-Wg = div{gF)-^{g*). 

□ 

The next result, Fedcrer's coarea formula, will be of critical importance to us in the sequel. 

11.2. Theorem ([Fed69], Theorem 3.1). Suppose X andY are Riemannian manifolds of dimension 
N and k respectively, with N > k. If f: X is a Lipschitz map, then 

(11.5) [ 9{y)Jf{y)dn''{y)= [ {[ g{y)dn^-'^{y)}dH''{x) 

whenever g: X ^M. is -integrable. Here, Jf{y) denotes the k-dimensional Jacobian of f at y, 

namely, the norm of the differential of f at y, df{y). Alternatively, it is the square root of the sum 
of the squares of the determinants of the k x k minors of the differential of f at y. 

11.3. Lemma. Let u: ^ M. be a Lipschitz function and let A C be a set of measure zero. 
Then 

W^"^(u"^(s) nA)=0 for almost all s. 
This can be directly seen from the coarea formula: 

= / \Vu{y)\dy= [ n'^-\Ar\u-\s))ds. 

One of the fundamental results of geometric measure theory is that any set of finite perimeter 
possesses a measure-theoretic interior normal which is suitably general to ensure the validity of the 
Gauss-Green theorem. 
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11.4. Theorem (DeGiorgi-Federer [DG61a, DG61b, Fed45, Fed58]). If E has finite perimeter, then, 
with V denoting the interior unit normal, 

f divFdy = - I F{y) ■ u{y)dn''-\y) 

JE Jd'E 

whenever F : is Lipschitz. 

The DeGiorgi-Federer result shows that integration by parts holds on a very large and rich family 
of sets, but only for fields, F, that are Lipschitz. As we explained in the introduction, the Gauss- 
Green formula for BV vector fields over sets of finite perimeter was treated by Maz'ja [BM69] and 
Volpert [Vol67a]. We contrast their work with that of the following result by Fuglede. 

12.1. Theorem (Fuglede [Fug60]). Let F e VMp{R^), 1 < p < do. Then there exists a function 

g : R-"^' R with g e L'^, ^ + ^ = 1, such that 

(12.1) / divF = - [ F{y) ■ u{y) dH^'-^y) 

JE Jd*E 

for all sets of finite perimeter, E, except possibly those for which 

[ g{y)dn''-\y) = ^. 

Jd'E 

The following, which is a direct consequence of Fuglede's result, will be of great use to us. 
Suppose that u: — > R is Lipschitz. For s <t, consider the "annulus" As;t := {x : s < u{s) < t} 
determined by u. Then, by appealing to the coarea formula, we see that As;t is a set finite perimeter 
for almost all s < t. Moreover, again appealing to the coarea formula, we see that, for almost all 
s <t, 

(12.2) {divF){A,,t):= I divF = - f F{y) ■ u{y) dn''-\y). 

One of the main objectives of this paper is to demonstrate that, when F G VM.'^iRJ^), we can 
extend Fuglede's result by showing that (12.1) and (12.2) hold for all sets of finite perimeter, not 
merely for "almost all" sets in the sense of Fuglede [Fug60] . Although we don't employ Fuglede's 
theorem directly, his result provided the motivation and insight for the development of our method. 

The Gauss-Green formula for bounded divergence-measure fields over sets of finite perimeter was 
first obtained in Chen- Torres [CT05]. The product rule from Lemma 11.1 was used to prove that 

(12.3) {divF){E^):= f divF = - [ 2xeF -Vue, 

Je^ Jd'E 

where 2xeF ■ Wue is the weak* limit of the measures 2xeF • Vuk- One of the main objectives of 
this paper is to obtain the trace measure as the limit of normal traces over smooth boundaries that 
approximate d*E. 

3. The Normal Trace and the Gauss-Green Formula for T>M°° fields over 

SMOOTHLY bounded SETS 

In this section we develop a method of obtaining the normal trace and thereby obtain the Gauss- 
Green formula for a bounded divergence-measure field over any smoothly bounded set. This provides 
the foundation for the development for the general case. 

12.2. Definition. Given a compact -manifold, M, we define the exterior determined by M to 
be that (connected) component, U, ofM.^ \M which is unbounded. The interior determined by M, 
U , is defined to be everything else in the complement of M ; namely, 

oo 

U = [j Bf;, Bk c R^ \ M a bounded component. 
fe=i 
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Thus, 

OO 

R'^ \M = UU { U Bk) = UUU. 

k=l 

13.1. Theorem. Let U C M.^ he the interior determined by a compact, , manifold M of dimension 
N—1 with Tl^~^{M) < oo. Then, for any F G VAi'^^{R^), there exist a signed measure a supported 
on dU = M such that a « 'H^~^ \_ dU and a function u : dU — > M such that 

(13.1) KU) ■■= {divF){U) = -a{dU) = - / (?, • i^){y)dn''-Hy) 

Jau 

and 

<C||F|| , 

M'^' Moo — II llcxj' 

where C is a constant depending only on N and U. 

With i/(y) denoting the interior unit normal to M at y, we may regard J'j • i/ as the interior 
normal trace of F on dU and thus write 

{3^,-i^){y) = F{y)-u{y). 

Hence, with this convention, it is convenient to abuse the notation and thus write (13.1) as 

(13.2) i^{U) = [ divF = - [ F{y) • i.{y) dH^'-^y), 

Ju JdU 

while bearing in mind that, since F is merely a measurable field and thus defined only up to a 
Lebcsguc null set, it may not even be defined on dU. We use the term "interior normal trace" to 
suggest that • is determined by the behavior of F in the interior determined by the manifold M. 
The proof will reveal that, in a similar way, it is possible to define the concept of "exterior normal 
trace". This will be discussed more fully below in Theorem 21.1. 

The next example shows that our general trace theorem remains valid even though the field F is 
singular on an open subset of dU. 

13.2. Lemma. Let F G T>A4'^^{'R^) whose distributional divergence is a measure /i, and let F^ be 
a mollification of F . Then, because Fg is smooth (in particular, Lipschitz), the classical divergence 
theorem holds whenever E g M.^ is a set of finite perimeter, namely, 

(13.3) / divF, = - / FM ■ u{y)dn''-'{y). 

JE Jd'E 

If, in addition, we assume the following two conditions: 

(i) Fe^F W^-i-a.e. ond*E, 

(ii) ^x{^E) = 0, 

then 

(13.4) n{E) = - [ F{y) ■ uiy) dH^'-'iy). 

The importance of this result is that, with assiunptions (i) and (ii), we obtain the Gauss-Green 
theorem for all sets of finite perimeter whenever is a bounded, measurable vector field with 
divF = fi. As stated earlier, our main objective is to obtain the same result without assuming (i) 
and (ii), by defining a suitable notion of normal trace for F on d*E. 

Proof. Since := divF^ ^ divF = in M{R^) and using the fact that Fg \_ d*E is uniformly 
bounded, we obtain from (13.3) that 

t,,{E) := / divF, = - [ FM . u{y) d'H^'-^y) ^ - / F{y) ■ 1/(2/) dH^'-^y), 

JE Jd'E Jd'E 

Hs{E) — > IJ-iE) (by assumption (ii)). 
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This establishes our result. □ 

Proof of Theorem 13.1. According to [Whi57], there exist a unit vector field V defined on dU 
and a number 6 > such that V is close to the interior normal on dU and that, with the line 
segment joining p G dU and q G U defined as 

A(p) ■.= {q:q = tVip) + (1 - t)p, 0<t<l}, 

then 

A*{p) ■.= X{p)nB{p,S)€U. 
We think of A* (p) as a quasi-normal to dU at p. Moreover, as p ranges over dU, the quasi- normals, 
A*(p), fill out a neighborhood U* of dU in a one-to-one way. That is, 

U* = [j A* (p) with A* (p) n A* (q) = when p q} 

The mapping n: U* ^ dU, which can be considered as the projection of U* onto dU along the 
quasi- normal. A*, defined by 

7r(gr) =p if g e A*, 
is clearly of class and thus, so is the mapping ?/* — > M defined by 

V'(9) •= k ~P\ where 7r(g) = p. 
Consequently, Sard's thc'orc'm implies that the critical values, V, of tjj arc of measure zero. Thus, the 
implicit function theorem implies that dUt := tp~^{t) is a manifold for almost all t G (0,(5) \ V. 
The manifold dUt can be considered as a deformation of dU along the line segments \{p). Thus, 
the sets Ut := {ip > t} arc open subsets of U with smooth boundaries and dUt = ip~^{t) for a.e. t. 
Observe that tp~^{'t) = {q & U : q = t{V{p) — p), p & dU). Since V € and dU is compact, it 
follows that V is Lipschitz with Lipschitz constant, say Cc, and therefore that 

(14.1) n^-\^-\t)) < Ccn^-\du). 

To see this, consider a local coordinate system on dU expressed as the inverse of its projection 
onto the tangent plane at a point p G dU. Consider Tp{dU) as a subspace in M.^ with the notation 
p = (p',0) e as identified with p' G W'^-K Thus, g: Tp{dU) dU will render dU as the graph 
of g so that g{p') = {p',y{p')), where y G and is the "height" function. Now define h: dU — > dUt 
by 

h{p',y{p')) = tV{p',y{p')) + (1 - tKp',y{p'j). 
Clearly, all the derivatives of ho g are uniformly bounded on compact subsets of Tp{dU) and then 

(14.1) is evident. An application of Lemma 13.2 shows that 

^Ji{Ut) = - [ F{y) ■ u{y) dH^'-^y) for a.e. t > 0, 

JdUt 

where /U := divF is a Radon measure. Thus, the function F ■ i/ \_ dUt induces a measure on R^. 
That is, for any Borel set B c K^, define 

at{B):= [ F{y).u{y)d'H''-\y). 

JBndUt 

Since F G and H'^-'^{dU) < oo, we see from (14.1) that n^-^{dUt) < C for some C > 0, which 
yields that the measures at, t > 0, form a bounded set in A1(R''^). Hence, there exist a sequence 
{tk} ^ and Radon measures a, cr+, and a~ with a = — a~ such that 

(14.2) (a+,<7,-,,atJ - (a+,a-,a) in A^(M^). 

Now we show that cr+ is supported on dU. On the contrary, let x G spt(cr+) \ dU and choose 
B{x,r) such that B{x,r) (1 dU = 0. Since x G spt(CT+), there exists (fi G C{B{x,r)) such that 
/ (pda~^ := cr+(<p) ^ 0. Then, since (fi is continuous, we find that crf{ip) — > ct+((^) ^ 0. This implies 



Here, we could have defined E* := (IR'^ \ dU) (1 E* to develop the notion of the exterior normal 
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that (Jti'-p) 7^ for all small i > 0, which leads to a contradiction since dUt fl B{x,r) = and 

Spt(<7+) C dUt. 

Clearly, 

(15.1) al{dUt,)^a+{dU), 
since 



whereas 



\immia+{dUt,) = liminf > = CT+{dU), 

fe— »(x> k—>oo 



linisupcr+(aC/tJ = limsupCT+(i^) < a+{K) = a+{dU) 



for any closed set K that contains C/, especially when K = M^. In the same way, we prove that 

(15.2) ai^{dUt,)^a-{dU), 
and hence 

(15.3) at,{dUt,) ^ aidU). 

Also note that the measure a is well-defined, which can be seen as follows: Let Ut'^ be another 
sequence of open sets with smooth boundaries defined by Uf^ := {ip > t^} for which Lemma 13.2 
applies. Then, assuming that tk > t'^. for all k, we have 

KUt',\Ut,)= [ divF-f divF 

= ~ I F{y) ■ u{y) dn''-\y) + / F(y) ■ i/(y) <m''-\y) 

Since Uf^ \ Ut^ is a monotone sequence of sets with Uf^ \ Ut^ 0, it follows that /u(C^f^ \ Ut^) — > 
and therefore that i^^t'^ ) ~ {dUt^ ) 0, which shows that a is well-defined. 

To show that a « W^-^ L dU, let A c dU with n^-'^{A) = 0. Then there exists an open set 

N 

G<sM^,GdA, such that G := \J Bk{rk) with Vrf"^ < s and n^-'^{G n dU) < e. Therefore, 
we have 



(G) < / \F{y) ■ «.(y)| dH^-i(y) < \\F\\^^^^^ H^-^G n M,) < eC ||F||^„ 
JGnM. 



'(G)' 



that is, for some C > 0, ||crt|| (G) < eC ||i^||^oo((3), which leads to the conclusion \\a\\ {A) = as 
desired. 

We note that, since the sets Ut^. are increasing, we have 

(15.4) ^(C/,J ^ ^(C/). 

Thus, using (15.3) and (15.4) and sending fc — »• oo in n{Ut^) = —at^ipUt^) yield 

(15.5) ii{U) = -a{dU). □ 

15.1. Remark. From the proof above, it can he seen that, if U were of class C^, then the interior 
normals to dU themselves would not intersect in a sufficiently small neighborhood of dU. 
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4. Almost One-Sided Smooth Approximation of Sets of Finite Perimeter 

We now proceed to establish a fundamental approximation theorem for a set of finite perimeter 
by a family of sets with smooth boundary essentially from the measure-theoretic interior of the set 
with respect to any Radon measure that is absolutely continuous with respect to Ti^^^. That is, 
we prove that, for any Radon measure /U on such that ji « H^~^, the superlevel sets of the 
mollifications of the characteristic functions of sets of finite perimeter provide an approximation by 
smooth sets which are ||/i||-almost contained in the measure-theoretic interior of E. This allows us 
to employ Theorem 13.1 which, after passage to a limit, leads to our main result. Theorem 21.1. 

16.0. Lemma. Let n he a Radon measure on M.^ such that fi << H^~^. Let E be a set of finite 
perimeter, and let Uk he the mollification ofX^- Then, for any t e (0, 1) and A^-t '■= {v '■ Uk{y) > t}, 

there exist e — e{t) and k* = k*{e,t) such that 



(i) 


M{Ak;t\E)<e 


for all < f < i and all k > k* . 


(ii) 


M{Ak.,t\E')<e 


for &\\\<t<l and all k>k* 


(iii) 


M{E'\Ak.,t)<e 


for all 5 < t < 1 and all k>k* 


(iv) 


y\\{E\Ak;t)<£ 


for all < t < i and all k > k* . 



Proof. We first show (ii). With i < t < 1, choose < e < t — ^. Let uj, denote (as usual) the 
mollification of X^. We know that Uk{y) Uj^{y) for 7i^~^-a.e. y and therefore the same is true 
for as well. By Egorov's theorem, for any e > 0, there is an open set C/g such that ||/i||(Ue) < s 
and that \uk{y) — Uj^{y)\ < e for all y and for all k > k* = k*{e,t). On A^.f \Us, we have 

t < Uk{y). 

Since Uk{y) < (y) + e on \Us, we have 

^<t-e< u^{y) u^{y) = 1 ^ y g E\ 

This yields 

Ak-t\U,(lE^ =^ Ak.t\E^cU,. 

Since ||yu||(i7£) < e, our desired result (ii) follows. 

For the proof of (i), given <t < |, we choose < e <t and proceed as above. 

We next show (iv). With < t < ^. choose < £ < ^ — t. For all large k, we have 
\ukiy) - < £ for all y ^ U^. Thus, on E\ U^, 

^ - Uk{y) < u^{y) - Uk{y) < £ Uk{y) >^-e>t fov all y G E\Us, 

which implies E\Us C Ak-t- Therefore, E \ Au-j. C and thus ||/i||(-E \ A^-^t) < £■ 

For the proof of (iii), given ^ < t < 1, we choose < £ < 1 — t and proceed as above. □ 

16.1. CoroUciry. For each <t < ^ and e < mm{t, ^ — t}, there exists k* = k*{e,t) > such that 

(16.1) |l/i||(Afe;tA£;) < £ whenever k > k* . 

For each ^ <t < 1 and e < min{t — i, 1 — t}, there exists k* = k*{s, t) > such that 

(16.2) \\n\\{Ak-tAE^) < e whenever k > k* . 

16.2. Remark. In the previous result, we used open superlevel sets Ak-i '■= {y : Uk{y) > t}. However, 
we could have used closed superlevel sets A^-^t := {y : Uk{y) > t} to obtain the same result. We also 
note that, for an arbitrary Radon measure w, we have 

(16.3) oj(Ak-,t) - co{Ak-,t) = oj{dAk-,t) = 



GAUSS-GREEN THEOREM, SETS OF FINITE PERIMETER, AND BALANCE LAWS 



17 



for all but countably many t, for the reason that the family of sets {dAj^.^ : t G M} is pairwise disjoint 
and any Radon measure w can assign positive values to only a countable number of such a family. 

17.1. Corollary. For each t > 0, there exist s{t) and k* = k*{e,t) > such that 

(i) ||/i||(^fc;tAi;) = ||/i||(Afe;tA£^) < s for all but countably many t € (0, |) and for k > k* ; 

(ii) \\n\\{Ak-t^E^) = \\^l\\{Ak■t^E^) < e for all but countably many t e (i, 1) and for k > k* ; 

(iii) \\n\\{dAk-t^E) = \\/.4{u~^{t)AE) < e for almost all t G (0, i) and for k > k* ; 

(iv) \\ii\\{dAk,tAE^) = \\ii\\{u'i;\t)AE^) < e for almost alltG (i,l) and for k > k* . 

For the case t = ^, only (i) and (iii) in Lemma 15.2 remain valid. To see this, we first show 

17.2. Lemma. Let jjL be a Radon measure on M.^ such that (j, « 7i^~^. Let E be a set of finite 
perimeter, and let Uk be the mollification of xe- Then, for t=\ and e > 0, there exists k* = k*{s) 
such that 

(17.1) M{E'\A,..)<e and M{A,..\E)<e. 

Proof. Since Uk{y) UE{y) for Ti^^^-a.c. y, the dominated convergence theorem implies that 
Wfc — >■ ue in L^(M^, ||/u||). Thus, given any £ > 0, for k large enough, we have 

(17.2) [ \uE-Uk\dM\> j {uE-Uk)dM>{l-\)M{E'\A^.^.), 

ft; 2 

which implies 

(17.3) M\{E^\A^.^.)<e. 
In the same way, we compute 

(17.4) !>/ \uE-Uk\dM>{\-Q)M{A^.^.\E), 
which implies 

(17.5) M{Ak..^\E)<e. 

□ 

The following remark shows that, with t = \ and with ji = TL^^^ L d*E > 0, (ii) and (iv) in 
Lemma 15.2 do not hold. 

17.3. Remark. // we define E := {y : \y\ < 1), then ul^{\) (ZM.^\E for all k, and therefore 
it is clear that 

(17.6) n^-\{A^.i\E^)nd*E)=n^-\d*E)^0 as/c^oo. 
// we now define E := {y G : \y\ > 1}, then u^^{^) C E for all k and thus 

(17.7) n^-\{E\Ak.i)nd*E) = n^-\d*E)^0 asfc^oo. 

17.4. Lemma. There exists a number < C < oo such that, for all positive integers k and almost 
allte (0,1), 

(17.8) n''-Hu^\t))<c. 

Proof. From Corollary 16.1, it follows that, for almost all t £ (0,1), there exists a sequence of 
smoothly bounded sets A^-^t such that either X^^ ^ X^i ;ti-a.e. (if ^ < t < 1), or X^^ ^ /z-a.e. 

(if < t < i). Since = H^''^ L d*E « H^-'^, it follows that X^^^ everywhere except 

for a set of Lebesgue measure zero. Since these functions are integrable, we may consider these 



18 



GUI-QIANG CHEN 



MONICA TORRES 



WILLIAM P. ZIEMER 



functions as distributions and thus, by appealing to Lebesgue's dominated convergence theorem, we 
may conclude 

X„ in the sense of distributions, 

and therefore 

Vx^ VX^ in the sense of distributions. 

Since all the functions are in BV, we find that 

VX,^^^ ^ VX, inX(M^), 

and the limit Vx^. is independent of t G (0,1). Therefore, the uniform boundedness theorem for 
measures (Corollary 4.2; also Lemma 4.3) implies that these measures are uniformly bounded in 
A^(M^); that is, there exists < C < oo, independent of t, such that, for almost all t € (0, 1), 

sup||VX,^J|(M^)<C. 

k 

Since 

l|VX,^J|=W^-^L«,-^(i), 
our result follows. □ 

Wc now offer another proof of this result here for the purpose of broadening the context of our 
development. By using the theory of integral currents [Fed69], the result is immediate. Rather than 
actually introducing integral currents, we will introduce a small structure that reflects the argument 
from integral currents. For this, let V denote the Banach space C^{K) of vector flelds ijj on K 
endowed with norm 

N 

m := sup(|V(2/)|+^|VVi(y)|), 

where K is a compact set such that E K. 
Let 



J E 

and, for almost every t e (0, 1), let 

Tk;tW:= i}{y)dy= i'{y)dy for V € V. 

J {y ■uk{y)>t} JAk-t 

Then, for ^ e V, we define the linear operators: 

dTsi^) := TE{divi>) = [ divi) dy = - [ i^-vdH^ 

Je Jd*E 



/N-1 



and 



dTk.4i^) := Tk.tidivi;) = [ div^ dy = - f i> ■ u dH''-\ 



where v is the interior unit normal. 
Since u'^^{t) is a C°°-manifold, then 

\\dTk;t\\--= sup \dn;tm = n''-\ui\t)). 



<1 



Indeed, with if) := defined on the manifold Mj, (t), the norm-preserving extension of tp to all of 

by Whitney's extension theorem yields the inequality 



\\dTk.,t\\:= sup \dTk.,tm>n''-\u-\t)). 

\\M\<1 



The opposite inequality is obvious. 
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Moreover, we find by the dominated convergence theorem that 

lim dn.t W ^ dT^ (V) for V € V, 

fc— ^oo 



and therefore, 



snp{\dTk,tW\} < oo for V G V. 
k 



By the uniform boundedncss principle (Theorem 4.1), we see that, since dTk-^t is a linear functional 
on V whose week limit, ST^, is independent of t, we have 

supH^ -\u-\t)) = sup WdTk-tW < oo, 

k k 

which gives our desired result. 

The above argument simply rephrases the following basic fact from the theory of currents. We 
know that, since E has finite perimeter, T^, is an integral current. Moreover, the currents Tj^-t 
converge to T^, weakly and therefore so do their boundaries, dTk-t — ^ 9Tj^; that is, 

for each smooth differential (A'' — l)-form a. Appealing to Corollary 4.2 yields our result. 

19.1. Lemma. Let u: f2 — > M he a Lipschitz function and E <s. Q a set of finite perimeter. Then 

n^-'^{d*E n = for almost all t. 

Proof. This result follows directly from Lemma 11.3, since we know that H^~^{d*E) < oo for any 
bounded set of finite perimeter, E c M^. □ 

19.2. Lemma. For almost every ^ < t < 1, we have 

(19.1) n^-\d*Enu^\t)) = o 

and 

(19.2) lim n^-^{d*E n Ak;t) = 0. 

Proof. This can be seen as follows. If we use Corollary 16.1 with fj, = H^~^ L d*E, we obtain 

lim fi{Ak-t \ E') = lim n''-^{Ak;t D d* E) = 0. 

k — >oo k — ^oo 

Clearly, (19.1) follows from Lemma 19.1 (see also Remark 16.2). □ 

19.3. Theorem (Approximation theorem). For almost every ^ <t <1, we have 

lim H^-\{E° U d*E) n u^\t)) = 0. 

k—^oo 

Proof. Since the Lebesgue measure is absolutely continuous with respect to 'H^~^, then using (16.2) 
in Corollary 16.1 with s> \ leads to 

|^fc;sA£;^| as A; -> 00. 

Therefore, since 

Rk;s := Ak;s \E^= Ak;s n (i;° U d*E°), 

it follows that 

(19.3) \Rk;s \ provided that s > ^. 
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Remark 16.2 indicates that we have the option of defining A^.t := {y : Uh{y) > t} without altering 
the development. With this option in force, we have u^^(s) C Ak^s and consequently, by the coarea 
formula (Theorem 11.2), 

jRk;3 Jo 

= [ n^-^{u^\t)n{E'^Lid*E°)) dt. 

J s 

Since Vuk VX^ and llVufeU^ < HVX^H (Lemma 7.1), it follows from Vitali's convergence 
theorem for s > 5 that 

/ IVwfel dy ^ 0. 

Thus, for a subsequence, we can conclude that, for a.e. t > s, 

n^-^{u^\t)n{E°Ud*E°)) ^0 ask ^00. 

The dependence on the subsequence is illusory. The reason is that, if there were a subsequence such 
that, for a.e. t, 

n^-'^{u^\t) n {E°Ud*E°)) ^a^O ask ^00, 
then one could appeal to our previous argument to conclude that, for some further subsequence and 
for a.e. t, 

n^-'^{u^\t)niE°Ud*E°)) ^0 asfc^oo, 
which is contrary to our assertion that a 7^ 0. 

Since s > 5 was fixed arbitrarily at the beginning of the proof, we conclude that, for a.e. t > |, 

W^-i {u^\t) n {E° U d*E°)) ^0 as A; ^ 00. 

□ 

5. Main Theorem 

In this section we establish our main result. Theorem 21.1. Let F G 'DA4'i^^{^Y). We define a 

measure ak;t for all Borel sets i? (s by 

(20.1) ak;t{B):= [ F{y) ■ u{y)dn''-\y), 

JBndAkft 

where F{y)-i/{y) denotes the normal trace of F on the smooth boundary dA^.f introduced in Lemma 
15.2. 

We begin with a lemma that will lead to several of the assertions in Theorem 21.1. 
20.1. Lemma. If E (S: ft is an arbitrary set of finite perimeter, then we have 

(20.2) [ F Vukdy = [ [ F -VkdH^'^dt 

Je Jo Jsnu^^it) 

for any F E L1^^{Q;R^), where Uk denotes the mollification of Xb introduced in Definition 6.3 
and Lemma 6.5. 

Proof. Let TV be the set on which Vu^ = 0. Then 

/ F-Vukdy= / F-S/ukdy+ F -Vukdy 
Je Jbxat Jm 

= [ \^Uk\^^^dy + 
Je\m 

= / \^Uk\gdy, 
Je 
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where g = ^'^'^ ^ ^^^ | — - . Then, by the Coarea Formula, we have 



[ F-Vukdy = I [ gdH 

Je Jo J u7^ (t)n(E\M) 



it)n(E\M) 



1 

N-l 



F-VkdW^-'dt 

Ju-\t)niE\^f) 
= [ [ F-VkdH^-^dt, 

Jo Ju7^(t)nE 



where we used Uk{y) = for y e u];\t) n {E \ M). □ 

With the help of Theorem 13.1 and the results in §4, we now establish our main theorem. 

21.1. Theorem (Main theorem). Let O C be an open set. Suppose that F e 'DM'^^(Q) with 
divF = n G A1(0). Let E <£ Q be a set of finite perimeter. Then 

(I) For almost every s e (^,1), there exist a signed measure Ui (independent of s) and a family 
of sets A^.g with smooth boundaries such that 

(a) M{Ak.sAE^)^0; 

(b) The measure Ui is the weak* limit of the measures ak-si 

(c) (7j is carried by d*E in the sense that \\(Ji\\ [Q. \ d*E) = 0; 

(d) \\ai\\ «n''-^Ld*E; 

(e) lim n^-\dAk:s n U d*E)) = 0; 

A;— *oo 

(f) lim \\ak.ME"ud*E)=0; 

k — ^oo 

(g) The density of Ui, denoted as ■ v, is called the interior normal trace relative to E 

of F on d*E and satisfies 



(y); 



(21.1) / divF=:ii{E^) = -ai{d*E) = -[ {"Ji ■ v){y) dU''-^ 

J E^ Jd*E 

(h) // {2F ■ Vuk)X^ is considered as a sequence of measures, then this sequence converges 
weak* to the measure {J'i ■ v)H^^^ \_d*E, i.e., 

{2F ■ Vuk)X^ ^ (J, • u)H^-^ L d*E in M{n)- 

(i) \\cJ^\\ = p^■u\\^.l,,E,n«--<\\noo■,E■ 
(II) For almost every s G (0, ^), there exist a signed measure (independent of s) and a family 
of sets Ak-s with smooth boundaries such that 

(a) MiAk-gAE) ^ 0; 

(b) The measure ag is the weak* limit of Uk-s', 

(c) CTe is carried by d*E in the sense that ||cre|| (O \ d*E) = 0; 

(d) llaell «H^-iL5*i5; 

(e) lim n^-^{dAk-s r\E)= lim n^-^{u7^{s) r\E) = 0; 

(f) lim ||afc;,||(0\£°)= lim \\cJk-A{E) = Q; 

k — >oo k — >QO 

(g) The density of a^, denoted as S'e • i^t is called the exterior normal trace relative to 

E of F on d*E and satisfies 

(21.2) / divF =: fi{E) = -ae{d*E) = - [ {J,- u){y) rfW^-^(y); 
Je Jd*E 

(h) // {2F ■ Wuk)X,^o is considered as a sequence of measures, then this sequence converges 
weak* to the measure [Te. ■ vyH^~^ \_ d*E, i.e., 

(2F • Vwfe)X^o ^ (S^e • i/)W^-i L d*E in M{Qy, 
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(i) llaell = llJ'e • t^\\oo;d'E,H«-^ < 1 1 1 1 oo;n\£; • 

Proof. Wc will prove only part (I), since the proof of part (II) is virtually identical. For notational 
simplicity, we will use the notation a rather than ai in the proof of part (I) . Throughout the proof, 
we will consider only those values of s € (^,1) for which the results in §4 are valid for all the mollified 
functions, u^, of thus omitting at most a set of measure zero. Without increasing the measure 
of this exceptional set, call it S, we will also include those values of s for which ||Mll(^fe ^(s)) 7^ for 
all k. For the rest of the proof, we fix such an s ^ <S. 

We start with (a). We consider the sets Ak-s as in Lemma 15.2. The desired result follows directly 
from Corollary 16.1. 

(b) Since F is bounded, Lemma 17.4 implies that there exists a constant C such that 

(22.1) hk;s\\m<C, 

which yields, as in (14.2), the existence of a signed measure o-g such that 

(22.2) ak-s ^ cTs in M{n). 

Utilizing (16.2), we also obtain that iJ,{Ak-s) Since Theorem 13.1 yields iJL{Ak-s) = 

— (Tfe;s(0), we obtain, after letting A; — > oo, 

(22.3) ix{E^) = -as{n). 

Since the left side of equation (22.3) is independent of s, we show next that CTs is also independent 
of s (and independent of the sequence in the weak* convergence (22.2)). To see this, we fix any 
(j) G C'c(^) ^"^d note that, since F is a divergence-measure field, the product rule in Lemma 11.1 
implies that (f)F is also a divergence-measure field. Proceeding as above with (j)F instead of F, we 
obtain 

(22.4) / div{(j)F) = - I (j)das 

for any (p € C^{n). Therefore, for any two measures as and as' that are the limits in (22.2), we have 
that (j>das = J^^ (pda^' for any (p G C^(0) and thus we conclude that as = as'. 

(c) Let A C \ d*E be an arbitrary Borel set. Referring to (5.4), we see that 

IIVXJI {A)=0. 

On the other hand, we know 

0=||VXJ|(A) 

(22.5) = inf{||VX^||(i7) : A C U, U open} 

= inf{||VxJ|(C/) ■.ACU,U open, ||VxJ| (OU) - 0}. 

In order to prove that ||a|| {A) = 0, we proceed by contradiction via assuming ||cr|| {A) > 0. From 

(22.5) , there is an open set U D A such that ||Vx^|| (dU) = and 

(22.6) IIVXJI {U) < tMl, 
Prom Lemma 20.1, we have 

(22.7) j \F-Vuk\dy= ( ( \F ■ Uk\ dU^-^ dt. 
Ju Jo Vi7n«~^(t) 
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Since U is open and ak-t ^ a in M{0,), 

\\a\\{A) < 2 / ||c7||(f/)di < 2 / liminf ||(Tfe.t||(C/)rft 

< 2/ liminf ||afe;t||(C/)(it <21iminf / \\ak;t\\ {U) dt, 
by Fatou's lemma. Therefore, we have 



\\a\\{A) < 21iminf/ / \F{y) ■ ,.iy)\ dH'^ -\y) dt 

Jo Ju-^{t)nu 

= 2Uminf / |F • Vufe| dy 

k^oo Jij 

< 2||F||^ lim / \Vuk\dy 
= 2||F||^||VxJ|(t/) 

< Ikll(^), 

where we used Lemma 6.5(iv) and the fact that ||VX^|| {dU) = 0. This yields a contradiction and 
thus establishes our result. 

(d) Let A C 1^ be a Borel set with = 0. Then, appeaHng to (5.4), we find that 
||VXe|| ip) = 0. From this, the proof can be proceeded precisely as in (c) to yield our desired 
conclusion. 

(e) This is the result of Theorem 19.3. 

(f ) In view of the definition 



ak-AB):= I F{y)-v{y)dn''-\y), 



and the fact that the normal trace is bounded, the result follows immediately from (e). 
(g) From (a) , we have the existence of smoothly bounded sets such that 

(23.1) ||Ai||(Afc;^A£;i) -> asfc^oo, 

where s > 1/2 was fixed at the beginning of the proof. From Theorem 13.1, we know that our 
desired result holds for smoothly bounded sets: 

(23.2) n{Ak;s):^ I divF = ~ [ F{y) ■ u{y) dn''-\y). 



We note that, with our notation in force, we may write (23.2) as 

(23.3) M(^fe;s) = -CTfe;s(0) = -ak-s{dAk-s)- 
Since 

(23.4) /i(Afe;s) — > /x(£;^) and o-fe;s(f2) ^ a-(f2) as fc ^ oo, 
we obtain 

/i(£i) = -a{d*E). 

Because ||c7|| << Ti^^^ L d*E, we know that there exists ■ u E L^{d*E) such that 



'Bnd'E 

which gives (21.1). 



<j{B)=l (Ji.u){y)dn''-\y), 

J Bnd'E 
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(h) From Lemma 20.1, we obtain 

lim / F-Vukdy= lim / / F{y) ■ Vk{y) dH^-'^{y) dt = lim / (Jk-t{E)dt. 



Thus 



X^F- Vu£;(0) := lim / X^F-Vukdy 



F ■ Vuk dy 

IE 



= lim / 

= lim / ak;t{E)dt+ lim [ ak;t{E°)dt (by (f) above) 
~ lim / afe;f(0)dt 

k^OG J j_ 

Let y be a function in C^(fi). Since (pF is also a bounded divergence- measure field, we can proceed 
as above with the vector field ipF instead of F to conclude that 

(24.1) [ ^dx^F-VuE = \ I ^da, 



which implies that a = iXj^F ■ Vue- 

(i) We have that, for W^-i-a.e. y G d*E, 

<^{B{y,r)) 



(24.2) "Ji ■ u{y) = lim 



r^o ||Vxs||(B(y,r))' 

where we can choose the balls B{y,r) such that ||Vx£;|| {dB{y,r)) = \\a\\ {dB{y,r)) = 0. 
Using a similar argument as in (c), we obtain 

\\aUBiy,r)) = 2 [\\aUB{y,r))dt 

= 2 lim / \\ak-t\\{B{y,r))dt 

fc— >CX) Jl 

= 2 lim [\\\ak■,t\\{B{y,r)^E^) + \\ak■,t\\{B{y,r)(^{E''i^^*E)))dt 

fe— >oo / 1 

2 

= 2 lim / \\ak;t\\{B{y,r)r^E^)dt 



2 lim / / \F-u\dn^-^dt 



k^ooji ^„-i(t)ns(j/,r)nBi 
<2j|F||^^^, lim / dH^'-'dt, 

where we have used the fact that \\crk;t\\ {E^ U d*E) ^ as fc ^ oo for a.e. t > 1/2. 
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Therefore, from (24.2), we obtain 

MiB{y,r)) 



\3'i-u{y)\ < lim 



IIVxeII iB{y,r)) 
< 2 IIFII^.,. lim lim 4 /^.-(^)nB(.,g '^^"-^^^ 

= 2||F|| lim lim K mnBiy,r) 

II ^^"^'^ r^ok^^ r _ .dH^-^dt 

JO Jui_ {t)nB(x,r) 

= ll-^lloo;Si = ll-'^llooiE- ^ 

As a direct result, we obtain the Gauss-Green theorem for divergence-measure fields over sets of 
finite perimeter. This also shows that our definition of the normal trace is in agreement with that 
given in the sense of distributions, Definition 7.3. 

25.1. Theorem (Gauss-Green theorem). Let Q C be an open set. Let F G 'DMf^^{0,) and let 
E <S:Q, be a bounded set of finite perimeter. Then, 

N-l 



(25.1) / ipdivF+ F-Vip^- ^p{^,-u)dn 

for all (fi G C^{n), where v is the measure-theoretic interior unit normal to E on d*E. 
Proof. From Lemma 11.1, it follows that (pF is a bounded divergence- measure field and 

(25.2) div{(fF) = pclivF + F -Vifi. 

Following the proof of Theorem 21.1 applied to ipF (instead of F), (see Theorem 11.1), we obtain 

(25.3) / div{ipF) = - [ {Ji ■ v) dn^-\ 

J E'^ Jd'E 

which, due to (25.2), gives the desired result. □ 

We conclude this section with the following remark. 

25.2. Remark. Theorem 25.1 implies that, when E is an open set of finite perimeter, our trace S^i-u 
agrees with the one defined in (8.1). 

6. The Divergence-Measures of Jump Sets via the Normal Traces 

In Theorem 21.1, we have defined the interior and exterior normal traces of G VJ\/['^^{Vl), 'Ji ■ v 
and 3^e • f , over the boundary of a set of finite perimeter E <sL In order to obtain the interior 
normal trace of F on where E := E^ \J d"''E, we reproduce the proof of Theorem 21.1 and 

apply it to E. Therefore, the trace measure, denoted by cr_, is obtained by using the level sets 
Bk;s = {vk > s} for some s G (5,1), where Vk is the mollification of Xg- We note that, for all y G fi, 

Pe * X^{y) +Pe* XM = 1> 

and therefore 

where 1 — s G (0, ^). Since —u is the interior unit normal to E, we have 

(25.4) a_{n) = - lim / F ■ u dU^'^ = - \mi ( F-vdU^''^ 

= - lim CTfe;i_s(M^) = -ae(0). 
The following observation now becomes evident. 
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25.3. Corollciry. The interior trace of F relative to E on d*E is the same as minus the exterior 
trace of F relative to E on d*E. 

In order to establish the relation between it, and cje, we subtract (21.2) from (21.1) and obtain 
the following formula for fi = divF. 



3.1. Corollary. n{d*E)= [ {Ji ■ - ■ i^){y) dH^'-^y). 

Jd*E 



26. 

'E 

We offer the following simple example to illustrate our result. This example also dramatically 
demonstrates the difference between the classical derivative and the weak (distributional) derivative. 

26.2. Example. Consider the most elementary situation: N = 1, Q. := (—1,2), E := [0, 1], and / is 

a non-dccrcasing function defined on (—1,2) which is continuous everjrwhere except at y = 0, 1, at 

which points wc assume that / has right-continuity. 

(i) Case i < s < 1. Since / is in BV , we know that f = n for some measure ji. Then, according 
to Theorem 21.1, 

li{E^) = M((0, 1)) := / f' = -Ji- - ■ 1.(0), 

where ■ v>{l) — limy_>i_ f{y) and S'i • i>'(0) — /(0+). Indeed, the sets Ak-s, with fixed ^ < s < 1, 
form a nested family of open intervals contained in [0, 1]. The measures cxfe correspond to / evaluated 
on the point masses located at yk'-, thus, as in (23.4), f{yk) converges to a limit, 3"; • v{l). 

(ii) Case < s < 5. Then the sets Ak-s, < s < 5, form a nested family of open intervals 
containing [0, 1]. Similar to the above, we have 

M(£;) = /x([0, 1]) := f = J,. - 1, ■ 1.(0), 
Jo- 

the measiu-os Uk correspond to / evaluated on the point masses located at yk, and thus "J^ ■ i'(l) = 
/(1+) and ■ i^(0) = lim^^o- f{y)- 

7. Consistency of the Normal Traces with the Classical Traces 

We now proceed to show the consistency of our normal traces with the classical traces when F is 
continuous. First we have the following lemma. 

26.3. Lemma. Let /j, = divF for F G VM°°{R^) n C{R^ ■,R^). Then 

M (G) = 

for any set G that can be written as the graph of a Lipschitz function f. 
Proof. First we have 

G:={{y',f{y')) : y'Gl^cK^-i}. 

By regularity of /x, it suffices to show that fi{K) ~ for any compact set K C G. Given any 
compact set K C G, let Uk C M.^ be a sequence of open sets satisfying 

(26.1) fi{UknG)^ fi{K). 

Fix any set Uk- We note by Besicovitch's theorem that Uk can be written up to a set of -measure 
zero as a countable union of disjoint open parallelepipeds I^ (the fact that we can use parallelepipeds 
instead of balls follows from Morse [Mor47]). Thus, we have 

(26.2) [jliCUk and M\ {Uk \ [j 1^ = 0. 

i=l i=l 

Denote Uk simply as U and as Ii. We fix an i and note that, for t small enough, the graphs 
Tt := {{y', f{y') +t) : y' &W ^ K^-i} and Bt := {{y', f{y') -t) : y' &W d M^"!} are contained 
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in /j. Let Rt be the region inside Ij, bounded above and below by and Bf respectively. For a.e. 
t, we define 



/ F{y)-u{y)dn''-\y) for a.e. 

J dB.A(TtUBA 



at 

JdBtMTtUBt) 

where i/{y) is the interior unit normal to Rt on dRt \ {Tt U Bt). Since Lemma 13.2 applies to Rt for 
a.e. t, we arrive at 

= / divF 

= - [ F{y) ■ 1/(2/) dn^'-Hy) - I F{y) • v{y) dU^'-^y) - a* 

JBt JTt 

= [ F{y',yr,-2t)-v{y)dn^-\y)- I F{y' ,y^) ■ u{y) dn^-\y) - at 

JTt JTt 

= f {F{y',yn-2t)-F{y',yn))-v{y)dn''-\y)-at. 

JTt 

Since F is continuous and at ^ as i ^ 0, wc find that there exists t^ie, F,G) > such that 

lx{Rt) < £ for all t < to{e, F, G). 



Then we have 



n G) = lim fj,{Rt) < £, 



which implies iJ.{Ii n G) = since e is arbitrary. Therefore, using (26.2), we obtain 
l^{K) = lim n{Uk n G) = lim V n G) = 0. 

□ 

27.1. Theorem. If F G VAAf^^ifl) is continuous and E Vt is a set of finite perimeter, then 
(Ji = F ■ VuE, and the normal trace "Ji ■ u is in fact the classical dot product F ■ u, where u is the 
interior unit normal to E. 

Proof We recall that, by definition, E = E^ d*E. Denote E = E^ d*E. Then we have 
F ■ VuE = lim / F-Vukdy 

/ xeF ■Wukdy+ ]iTD. / Xs^-'^Ukdy. 

' JRN JRW 



= lim 

k — ^oo 



If Vk denotes the convolution Xe * Pi/k, since Wfe + t;fc = 1, we obtain 

F ■ VuE = lim / XeF ■ Vuu dy - lim / Xe^ ' ^'^k dy 

= y + y = i(-. + <^i-Ma*i;)) 

= - ^fi{d*E), 

where we used Theorem 21.1 (h) and Corollary 26.1. 

Since d*E is an {N — l)-rectifiable set (see (5.5)), it follows from Lemma 26.3 that 

M {d*E) = 0, 

that is, 

F ■ Vue = cTj. 
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Thus, for H^-i-a.e. y G d*E, 

r^o \\\7xE\\{Biy,r)) r^Ok^^ J^^^^^^d\\VxE\\ 

Since Vuk — > Vxb weak* and F is continuous, and noting that rj can be chosen such that 
||Vxb|| {dB{y,rj)) = 0, we obtain 

(J, • = hm r = hm — ^^^^^ — — -— 

= F{y)-u{y), 

by differentiation of measures. □ 

The following corollary gives more information of the trace cTj and the level sets u'^^{s) when 
1 

2' 

28.1. Corollary. The trace measure di given in Theorem 21.1 satisfies 

ai{U.^) = 2 lim / F • Vufe = 2 lim lim / F-Vukdy 



S i + - 



= 2 lim / F-Wukdy. 

k^ooJA, 1 



Proof. Theorem 21.1 (h) shows 

(28.1) ai{R^) = 2 lim / F Vuk dy. 
Using Lemma 20.1, we find 

c7i(M^) = 2 lim / CTi(M^)di = 2 lim lim / / F-vdH^''^ 

= 2 hm lim / F ■ Vuk dy. 



One can easily verify that the limits s ^ \+ and k ^ oo can be interchanged. Noting that 
(J Ak-s = Af.. 1 , we conclude 

cri(IR^) = 2 hm hm / F Vwfe dy = 2 lim F Vuk dy. 

□ 

8. One-sided approximation of sets of finite perimeter 

It is well-known that, a set of finite perimeter, E, cannot be approximated by smooth sets that 
lie completely in the interior of E. For example, consider the open unit disk with a single radius 
removed, and let U be the resulting open set. Then the Hausdorff measure of the boundary of U is 
27r plus the measure of the radius, while the Hausdorff measure of the reduced boundary is 27r. Thus, 
if Uk is an approximating open subset of U, then its boundary will be close to that of boundary U 
and so its the Hausdorff measure will be close to 2?: plus 1. Adding more radii, say m of them, will 
force the approximating set to have boundaries whose Hausdorff measure close to 2^ plus m. In 
general, if we let K denote any compact subset without interior and of infinite Hausdorff measure, 
then the approximating sets will have boundaries whose measures will necessarily tend to infinity. 
On the other hand, the one-sided approximation is possible for open sets of class , see Theorem 
13.1. More generally we have the following: 
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28.2. Proposition. Let U C M.^ be an open set with Tl^~^{dU) < oo. Then there exists a sequence 
of bounded open sets Uk C Uk C U such that 

(i) \Uk\ = m,- 

(ii) \Uk\ - \U\; 

(iii) n^-^idUk) -^n^-Hdu). 

Proof. By definition, for eacii integer k, there exists a covering of dU by balls 
each with radius r^, such that 

oo oo ^ 

i=i i=i 

where u)n-i is the H^~^ measure of the boundary of the unit ball in M^. Since dSl is compact, the 
covering may be taken as a finite covering, say by m of them, Bi{ri), Bi{r2), • . • , B^{rm)- Then 
the open set Vk := [j Bi{ri) has the property that 

m 

dVk c U dBiin) 

i=l 

and therefore that 

oo 

m 1 

'W'-'iOVk) < U dBi{n)) < Va;^-irf-i < H^^-^dU) + ^. 

Thus, the open sets Uk ■= U \Vk C U will satisfy our desired result, except that they are not 
smooth. □ 

Given an arbitrary set of finite perimeter, E, we know from §4 that E can be approximated by 
sets with smooth boundaries essentially from the measure-theoretic interior of E, that is, a one-sided 
approximation can "almost" be achieved (see Theorem 21.1(e)). On the other hand, the next result 
shows that, if E is sufficiently regular, there does, in fact, exist a one-sided approximation. The 
condition of regularity we impose is similar to Lewis's uniformly fiat condition in potential theory 
[Lew88a]. 

29.1. Theorem. Suppose that E is a bounded set of finite perimeter with the property that, for all 
y G dE, there are positive constants cq and ro such that 

\E°nB(y,r)\ 

29.1 ' ,1 ^' > CO for all r < tq. 

\B{y,r)\ 

Then there exists t G (0, 1) such that 

(29.2) Ak;t € E for large k. 

Proof. Choose a mollifying kernel p such that p = 1 on B{0, i). If 2/ € dE, we have 

\ (' X — y 

vkiy) ■■=X^N^^* Pedv) = — / \«XEix)p{— — )dx 

^ -TV / _ >^R",^(a;)da; 



m^\E)nB{y,f)\ 



\E^nB{y,^)\ ^ 

-> Co/2" := Co, 
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where < Cq < 1 depends only on the dimension N and is independent of the point y. Note that 



Thus, taking 1 — cq <t< 1, we see that Ak-t H dE = 0. Consequently, each connected component of 
the open set Ak-t lies either in the interior of E or in its exterior, and thus must lie in its interior. □ 

30.1. Corollary. Let E be a bounded set of finite perimeter with uniform Lipschitz boundary. Then 
there exists T G (0, 1) such that A^-t <e E. 

Proof. Since E has a uniform Lipschitz boundary, for each x G dE, there is a finite cone, Cx, with 
vertex x that completely lies in the complement of E. Each cone Cx is assumed to be congruent 
to a fixed cone C. This implies that the hypothesis of Theorem 29.1 is satisfied. Therefore, there 
exists < T < 1 such that Uk{y) < T for all k and all y G dE. □ 

30.2. Definition. An open set U C is called an extension domain /or € 'DAd°°{U), if there 
exists a field F* € 2?A^°°(R^) such that F = F* on U. 

30.3. Theorem. An open set U satisfying T-L^~^{dU) < oo is an extension domain for any F € 
'DJv[°°(U). More generally, if F & 'DAi'^^{^), then any open set of finite perimeter U Cl is an 
extension domain for F. 

Proof. We define an extension of F by 



We consider first the case 'H^~^{U) < oo. Let Uk be the sequence of approximate sets given in 
Proposition 28.2. Therefore, for any (p G C^(]R^) with \(p\ < 1, we employ our general Gauss-Green 
theorem, Theorem 21.1, to obtain 



Uk{y) + Vk{y) = 1 for all 2/ G R^. Therefore, for all y G dE, 



Ukiv) = 1 - Vk{y) < 1 - Co. 



F*{y) — \,{y)F{y) for all y G K' 



According to Definition 7.3, it suffices to show that 





Thus 




< 



\\divF\\{Uu) + \\F\\^n''-\dUk) 



< 



\\diyF\\{U) + \\F\\^n''-\dUk). 



Letting fc ^ oo, we obtain 




Thus, 
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We now consider the case that U (g is a set of finite perimeter and F G I?A^^^(r2). Proceeding 
as above and using Theorem 21.1, 

/ F*-Vipdy= I F*-V^dy 

= - [ ipdivF- [ (^Ji-i^dW^-i 

Ju Jd*U 

< \\divF\\{U) + \\F\\^n^-'^{d*U) < 00. 

□ 

31.1. Corollary (Chen-Torres [CT05]). Let U CR'^ be a bounded, open set with H^-'^{dU) < oo. 
Let Fi e VM°°{U) and F2 € VM"'{R'^ \U). Then, with 



F{y) 

we have 



Fi{y) yeU, 

F2{y) 2/eM^\C/, 



F e VM°°{m.^). 

Proof. Applying the previous result to 

Fl := X„{y)Fi{y) for all 2/ e 

and 

Fl := \«,^F2(y) for all y G M^, 

we see that 

F = F\ + Fl □ 



9. Cauchy Fluxes and Divergence-Measure Fields 
The physical principle of balance law of the form 

(31.1) / f{y,u{y))dH''-\y)+ j b{y)dy = Q 

JdE Je 

is basic in all of classical physics. Here, v>{y) is the interior unit normal to the boundary dE of E. 
In mechanics, / represents the surface force per unit area on dE, while in thermodynamics / gives 
the heat flow per unit area across the boundary dE. 

In 1823, Cauchy [Cal823] (also see [Cal827]) established the stress theorem that is probably the 
most important result in continuum mechanics: If both f{y,i^{y)), defined for each y in an open 
region Q and every unit vector v, is continuous in y and b{y) is uniformly bounded on Q, and if 
(31.1) is satisfied for every smooth region E (s i}, then fiy^f) must be linear in v. The Cauchy 
postulate states that the density flux / through a surface depends on the surface solely through the 
normal at that point. For instance, if f{y, v) represents the heat flow, then the Stress theorem states 
that there exists a vector field F such that 

f{y,u) = F{y)-v. 

Since the time of Cauchy's stress result, [Cal823, Cal827], many efforts have been made to 
generalize his ideas and remove some of his hypotheses. The first results in this direction were 
obtained by Noll [Nol59] in 1959, who set up a basis for an axiomatic foundation for continuum 
thermodynamics. In particular, Noll [Nol59] showed that the Cauchy postulate may directly follow 
from the balance law. In [GM76] , Gurtin-Martins introduced the concept of Cauchy flux and removed 
the continuity assumption on /. In [Zie83], Ziemer proved Noll's theorem in the context of geometric 
measure theory, in which the Cauchy fluxes were first formulated at the level of generality with sets 
of finite perimeter in the absence of jump surfaces, "shock waves". 
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However, as we explain below, all the previous formulations of (31.1) do not allow the presence 
of "shock waves" ; one of our main intentions in this paper is to develop a theory that will allow the 
presence of "shock waves" . 

In this section we first introduce a class of Cauchy fluxes that allows the presence of the exceptional 
surfaces or "shock waves" and then prove that such a Cauchy flux induces a bounded divergence- 
measure (vector) field F so that the Cauchy flux over every oriented surface can be recovered through 
F and the normal to the oriented surface. Before introducing this framework, we need the following 
definitions. 

32.1. Definition. An oriented surface in fl is a pair {S,i/) so that S fl is a Borel set and 

V : §^^^ is a Borel measurable unit vector field that satisfy the following property: There is 

a set E <s fl of finite perimeter such that S C d*E and 



where xs is the characteristic function of the set S and (y) is the interior measure-theoretic unit 
normal to E at y. 

Two oriented surfaces {Sj,i'j),j = 1,2, are said to be compatible if there exists a set of finite 
perimeter E such that Sj Cd*E and J^jiij) = (y) Xs^- (y), j = 1,2. For simplicity, we will denote 
the pair [S, v) simply as S, with implicit understanding that S is oriented by the interior normal of 
some set E of finite perimeter. We define —S = {S, —if), which is regarded as a different oriented 
surface. 

32.2. Definition. Let fl be a bounded open set. A Cauchy flux is a functional T that assigns to 
each oriented surface S := {S, i/) (s O a real number and has the following properties: 

(i) T{Si U 5*2) = J^{Si) + ^(82) for any pair of com,patible disjoint surfaces 81,82 <e f^; 

(ii) There exists a nonnegative Radon m,easure a in such that 



for every oriented surface S* d f2 satisfying C7{8) = 0. 

This general framework for Cauchy fiuxes allows the presence of exceptional surfaces, "shock 
waves" , in the formulation of the axioms, on which the measure a has support. On these exceptional 
surfaces, the Cauchy flux has a discontinuity, i.e., J^{8) 7^ —T{—S). In fact, the exceptional 
surfaces are supported on the singular part of measure a. When a reduces to the TV-dimensional 
Lebesgue measure C'^ , the formulation reduces to Ziemer's formulation in [Zie83] and in this case a 
vanishes on any 7Y^~^-dimensional surface, which excludes shock waves. 

The theory developed in this paper allows to approximate the exceptional surfaces or "shock 
waves" with smooth boundaries and rigorously pass to the limit to recover the flux across the shock 
waves. This allows to capture measure production density in the formulation of the balance law 
and entropy dissipation for entropy solutions of hyperbolic conservation laws. Once we know the 
flux across every surface, we proceed to obtain a rigorous derivation of nonlinear systems of balance 
laws with measure source terms from the physical principle of balance law in §10. The framework 
also allows the recovery of Cauchy entropy fluxes through the Lax entropy inequality for entropy 
solutions of hyperbolic conservation laws by capturing entropy dissipation; see §11. 

The main theorem of this section is the following. 



32.3. Theorem. Let T be a Cauchy flux in Ct. Then there exists a unique divergence-measure field 
F e VM^^{n) such that 



t^{y) = ^Biv) xs{y), 



(in) 



\T{d*E}\ < a{E) 

for every set of finite perimeter E <S:fl satisfying a{dE) = 0; 
There exists a constant C such that 



\^{8)\ <cn^-\8) 



(32.1) 
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for every oriented surfaee (S*, u) d fl, where ■ u is the normal trace of F to the oriented surface. 

When a reduces to the A^-dimensional Lebesgue measure as in Ziemer's formulation, the 
vector field F satisfies divF e L°° and J^{S) = —J^{—S) for every surface S, which thus excludes 
shock waves where the Cauchy flux T has a discontinuity, i.e., T(S) ^ —!F{—S). 

In order to establish Theorem 32.3, we need Lemmas 33.2-35.1 that were first shown in Degiovanni- 
Marzocci-Musesti [DMM99]. Here we offer simplified proofs of these facts for completeness. In 
particular. Lemma 33.2 is in fact a direct application of Theorem 33.1 (due to Fuglede) below, and 
Lemma 35.1 follows by an approximation and Theorem 21.1. We also refer to Schuricht [Sch07] for 
a different approach in formulating the axioms in Definition 32.2. 

The following theorem, due to Fuglede, is a generalization of Riesz's theorem, whose proof can 
be found in [Fug55]. 

33.1. Theorem. Let be a nonnegative measure defined on a a -field V of subsets of a fixed set 
X and X e V. Let y> be an additive set function defined on a system of sets U C V such that all 

finite unions of disjoint sets from U, together with the empty set, form a field J- which generates V. 
Assume that IJ,{A) < oo for every A cU. Then there exists a function g{y) G L^{X,V,ijl) with the 
property that 



if and only if the following hold: 

(i) For every e > 0, there exists 5 > Q such that Y^=i — ^ Z^*" every finite system of 
disjoint sets Ax, A^, An from U for which Y^^=i l^i^i) < ^! 

(ii) There is a finite constant C such that Ifi^i)] ^ ^ for every finite system of disjoint sets 
Ai,A2, ...,An from U. 

The function g is then essentially uniquely determined. Under the additional assumption that l-i{X) < 

oo, condition (ii) is a consequence of condition (i). 

Let {/} be the collection of all closed cubes in M.^ of the form 



where ai, bi, a2, ^2, cln, and b^ are real numbers. For almost every tj e [aj,bj], we define 



We define the vectors ei, 62, ejv so that the j-th component of ej is —1 and the other components 
are zero. We orient the surface 1^. with the vector ej. 

33.2. Lemma. Let T be a Cauchy flux in ^l. Then there exists a divergence-measure field F € 
^■^ioc(^) ^^c/i that, for every cube I = [ai,bi] x [02,62] x ... x [aAr,6jv] <£ and almost every 




I = [ai,6i] X [02,62] X ...[ajv,6jv], 



hi ={y&i-yj= Tj}. 



Tj e [aj,bj]. 




Proof. Step 1. We fix j G {1, N}. For every cube 7 C O, we define 




We have 



(33.1) 




Thus, from Theorem 33.1, there exists a function f^ e 



(O) such that 




for every /. 
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In fact, inequality (33.1) implies that G L°°{Q) since, for £^-a..e. y, 

■' \i\^o,yei \I\ - 

Fubini's theorem implies that 

(34.1) /x^(J)= r{Ir,)dTj= f Pdy= r / PdU^'-^dTj. 

Let Tj G [aj,bj], akj, and pkj be sequences such that 

ak,j <Tj < Pk,j , 

where akj is an increasing sequence that converges to tj as fc ^ oo, and (3k,j is a decreasing sequence 
that converges to tj as A; — > oo. Thus, from (34.1), we obtain 

(34.2) \ P HIr,)drj = ^ P [ f dn^'-'dr,. 

ak,j - Pk,j J a, O^kJ - Pk,j Ja^ Jl^. 

We let fc ^ oo to obtain that, for a.e. r, , 

(34.3) J^{Ir,) = J fdU'^-K 

Define 

F.= {f\f,...,f^). 
Then we find that, for almost every tj, j & {1,2, ...,N}, 

(34.4) HIr,) = -f Fiy).ejdn''-\y). 

Step 2. We now prove that the divergence of F, in the sense of distributions, is a measure. We 
define, for a.e. cube I, 

(34.5) m--^-f F{y).u{y)dn''-\y), 

Jdl 

where v is the interior unit normal to dl. 

Prom Step 1 and the definition of Cauchy flux, we have 

(34.6) m\ = \HdI)\<<r{I) 

for almost all closed cubes. Thus, we can again apply Theorem 33.1 to conclude that there exists a 
function g £ L^{fl; a), uniquely defined in O up to a set of a-measure zero, such that 

(34.7) 77(7) = - / F{y) • u{y) dH^-'iy) = [ g{y) da 

Jdl J I 

for almost every closed cube 7 c fi. 

Denote a the measure given by gda in Q.. We now prove 

(34.8) div F = 5- 

in the sense of distributions in any open set C/ (1 f7. 

Let / (E [/ be any closed cube. Then, for any G with support contained in 7, 



dy. 



(34.9) / F ■y(j)dy =Viui / • V$!i = - lim / (/)divFe 

where F^ = F * and p is the standard mollifying kernel. We now prove that, for £^-a.e. y 

ae{y) = divF^(y), 
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where is the convolution of function with the measure a; that is, 

(34.10) dM := (pe * ^){V) = I Ps{y - x) da{x). 

Jq 

Prom (34.6)-(34.7), we find that, for s < dist{dU,dn), 

I divFe{y)dy = - F,{y) ■ u{y) dy = - / F{y - x) ■ u{y)p,{x) dxdy 
J I Jai Jdl JR" 

= - / F{y - x) ■v{y)pe(.x)dydx 

JR" Jdl 

= - ( F{y) ■u{y)dy)pe{x)dx 

= / HIx)psix)dx, 
Jr" 

where = {y ■ cii < yi — Xi < bi, i ~ I, ■ ■ ■ , N}. We can consider the smooth function p^ as a 
measure in R^, say A^, by defining \e{^) = JaP^^^)*^^ Borel set A. We can also extend 

the measure d by zero outside Vt. Therefore, we find 

(35.1) / a{I^)pe{x)dx = {a*X,){I) = {X,*a){I)= [ X,{I^)da{x). 

Prom (35.1) and using (34.10), we compute 

/ Xe{Ix)da{x)= / Xe{Ix)da{x) = ( pe{y)dy)da{x) = / pe{y - x) dyda{x) 
jR'v Ja jQ ^Ji^ ' Jo. J I 

= (^j^Ps{y-x)da{x)jdy = j {pe * d^){y) dy 
= j^^e{y)dy. 



Therefore, 



divFe{y) '^y = ^e(y) '^y- 



Since the cube / g [/ is arbitrary, this shows that CFe{y) = divF^iy) for £^-a.e. y & U. Using this 
in (34.9), we obtain 

(35.2) / F • V(/) = - lim / (pdivF^ dy = - lim / cpa^dy ^ - / (j){y) da{y), 

Ju ^^°Ju ^-*^Ju Ju 

since the sequence of measures converges locally weak* in to a as e — > 0. □ 

35.1. Lemma. Let T be a Cauchy flux in Q. Then there exists a unique divergence-measure field 
F G VM'i^^in) such that 



(35.3) J^{S) = - / Ji-vdn 



^(5) = - / Ji-ui 
Js 



jN-l 



for almost every oriented surface {S, u) g Q,; that is, every surface S in fl satisfying the condition 
that a{S) = 0. 

Proof. Using Lemma 33.2, it follows that there exists an e 'DMi^^{Q) such that, for any cube 

I = [ai, 6i] X [a2, 62] X ... x [gn, 6^] <s ^, 

(35.4) ^(/.^ ) = - y F(y) • e, dU^'-'iy) 

for almost every Tj € [aj, 
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Let {S, u) be an oriented surface satisfying <j{S) = 0. Then, since divF = gda and proceeding as 
in Lemma 26.3, we conclude that ||divJ^||(6') = 0. We approximate S with closed cubes such that 



(35.5) 



where each is a finite union of closed cubes and Jj+i C Ji- Since S is an oriented surface, there 
exists a set of finite perimeter, E, such that S C d*E. Using Lemma 6.2, we have 

(36.1) d*{Jir\E) = SUid*Jir\E)UAfi, 

where lim 'H^-'^{Ni) = and thus, since J^{J\fi) < CH^-^{J\fi), we obtain 

(36.2) lim T{J\fi) = 0. 

The definition of Cauchy flux implies that 

\J^{d*{JinE))\<a{JinE), 

and thus (35.5) implies that 

(36.3) lim \J^{d*{Ji nE))\< a{S) = 0. 

z— >-oo 

On the other hand, using Theorem 21.1, we have 



(36.4) 

which yields 
(36.5) 



d'iJiHE) 



lim 



divF 



JiDE 



< \\divF\\ {JiHE), 



J 



N-l 



< \\divF\\ {S) = 0. 



Using (36.1)-(36.3) and Lemma 33.2, we obtain 

(36.6) lim \J^{d*{JinE))\ = lim \J^{S) + J^{d* Ji n E)\ 



= lim 

i—^oo 



F{y)-u{y)dn 



N-l 



d'JiHE 



0. 



From (36.1), (36.2), and (36.5), we obtain 



(36.7) 



lim 



/ F{y) ■ u{y) dn"-' - [ li-u dH"-' 

Jd'JinE Js 



0. 



Combining (36.6) with (36.7) yields 



^i-vdH 



N-l 



Assume now that there exists another vector field G = {g^,g'^, ...,g^) such that (35.3) holds. 
Then, for fixed j G {1, 2, .... N}, we have 

(36.8) / rdy= I fdn''-'dTj= r / g^dn''-UTj= [ g^ dy 

for any cube /. This implies that 

fHy) = g''{y) ^or almost every y. 

□ 

With Lemmas 33.2-35.1, we now prove Theorem 32.3 to explain how the Cauchy flux can be 
recovered on the exceptional surfaces based on the theory established in §3-§8. 
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Proof of Theorem 32.3. Let {S, u) be an oriented surface on which <j{S) ^ 0, i.e. J^{—S) ^ —J^{S). 
By definition of oriented surfaces, there exists a bounded set of finite perimeter, E := E^\Jd*E g O, 
such that 

(37.0) S^d*E and v{y) = UE{y)xs{y), 
where i^siy) is the interior normal to at y G 5. Consider 

E = E°U d*E. 

Then Theorem 21.1 implies that there exist the normal traces 3^, • v and J'e " defined ond*E = d*E 
respectively such that 

/ divF = - [ Ji-vdH^-'^, 

[ divF = - [ ye-vdn^-^ = - ( %-iydn'^-\ 

Je Jd'E Jd'E 

Theorem 21.1 indicates that the traces Ji ■ v and 7^ ■ v can be recovered, up to a set of arbitrary 
small 7i^~^-measure, from the neighborhood behavior of the vector field F . This observation allows 
us to define 

(37.1) j:{S)=T{S,u):=- [ Ji-vdH'^-\ 

Js 

and 

(37.2) jr(-5) = jr(5', -I/) := - / 3^e-{-u)dH^-^ = [ Je-i^dn^-\ 

Js Js 

In this way, we can recover the Cauchy flux !F through the corresponding divergence-measure field 
F over all oriented surfaces, especially including the exceptional surfaces. That is, the normal traces 
of F e I?7W^g(n) are the Cauchy densities over all oriented surfaces. This completes the proof of 
Theorem 32.3. □ 



10. Mathematical formulation of the balance law and Derivation of Systems of 

Balance Laws 

In this section we first present the mathematical formulation of the physical principle of balance 
law (31.1). Then we apply the results established in §3-§9 to give a rigorous derivation of systems 
of balance laws with measure source terms. In particular, we give a derivation of hyperbolic systems 
of conservation laws (38.7). 

A balance law on an open subset n of M.^ postulates that the producfAon of a vector-valued 
"extensive" quantity in any bounded measurable subset (s with finite perimeter is balanced by 
the Cauchy fl,ux of this quantity through the measure-theoretic boundary d^E of E (see Dafermos 
[Daf98, Daf05]). 

Like the Cauchy flux, the production is introduced through a functional P, defined on any bounded 
measurable subset of finite perimeter, E dVL, taking value in M*^ and satisfying the conditions: 

(37.3) V{Ei\jE2)=V{Ei)+V{E2) ifEinE2 = $, 

(37.4) \P{E)\ < a{E). 

Then the physical principle of balance law can be mathematically formulated as 

(37.5) J^id^'E) = V{E) 

for any bounded measurable subset of finite perimeter, E cfl. 
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Fugele's theorem, Theorem 33.1, indicates that conditions (37.3)-(37.4) imphes that there is a 
production density P G j'W(0;R*^) such that 



(37.6) V{E) = [ P{y). 



On the other hand, combining Theorem 21.1 with the argument from §9, it follows that there 
exists F e VM'^^{Vt-R^'''^) such that 

(38.1) J^id^'E) = - [ (Ji • I/) dn^-'^ = [ div F{y) 

for any set of finite perimeter, E <^ ^l. 

Then (37.5)-(38.1) yields the following system of field equations 

(38.2) divF{y) = P{y) 

in the sense of measures on fl. 

We assume that the state of the medium is described by a state vector field u, taking value in 
an open subset U of M*^ , which determines both the flux density fleld F and the production density 
field P at the point y £flhy the constitutive equations: 

(38.3) F{y) := F{u{y),y), P{y) := P{u{y),y), 

where F{u, y) and P{u, y) are given smooth functions defined on ?7 x Q. 

Combining (38.2) with (38.3) leads to the quasilinear first-order system of partial differential 
equations 

(38.4) divF{u{y),y) = P{u{y),y), 

which is called a system of balance laws (cf. [Daf98]). 
liV = 0, the previous derivation yields 

(38.5) divF{u{y),y)=0, 

which is called a system of conservation laws. When the medium is homogeneous: 

F{u,y) = F{u), 

that is, F depends on y only through the state vector, then system (38.5) becomes 

(38.6) divF{u{y)) ^ 0. 

In particular, when the coordinate system y is described by the time variable t and the space 
variable x = {xi, - ■ ■ , a;„): 

y ^ (t,xi,- ■ ■ ,Xn) = it,x), N = n + 1, 
and the flux density is written as 

F{u) = (u, /i(m), • • • , /„(«)) = {u, /(m)), 
then we have the following standard form for the system of conservation laws: 

(38.7) + Vx ■ /(w) = 0, xeM", ueK'=. 
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11. Entropy Solutions of Hyperbolic Conservation Laws 

We now apply the results established in §3 - §9 to the recovery of Cauchy entropy fluxes through 
the Lax entropy inequality for entropy solutions of hyperbolic conservation laws by capturing entropy 
dissipation. We focus on system (38.7) which is assumed to be hyperbolic. 

38.1. Definition. A function : M*^ — > M is called an entropy of (38.7) if there exists g : — > M" 

such that 

(38.8) Vqj {u) = VTi{u)Vfj {u), j = 1, 2, . . . , n. 

Then the vector function q{u) is called an entropy flux associated with the entropy r](u), and the 
pair {ri{u),q{u)) is called an entropy pair. The entropy pair {r]{u),q{u)) is called a convex entropy 
pair on the domain U C M*^ z/ the Hessian matrix V^77(u) > for any u £ U. The entropy pair 
{r]{u),q{u)) is called a strictly convex entropy pair on the domain U ifV^r]{u) > for any u G U. 

Fricdrichs-Lax [FL71] observed that most of systems of conservation laws that result from contin- 
uum mechanics are endowed with a globally defined, strictly convex entropy. The available existence 
theories show that solutions of (38.7) are generally in the following class of entropy solutions. 

39.1. Definition. A vector function u = u{t,x) G Lf^^{R^ x M") is called an entropy solution if 
u{t, x) satisfies the Lax entropy inequality: 

(39.1) dtr]{u{t, x)) + • q{u{t, x)) < 

in the sense of distributions for any convex entropy pair [t], q) : R*^ ^ R x M". 

Clearly, an entropy solution is a weak solution by choosing r]{u) = ±u in (39.1). 

One of the main issues in conservation laws is to study the behavior of entropy solutions in this 
class to explore to the fullest extent possible all questions relating to large-time behavior, uniqueness, 
stability, structure, and traces of entropy solutions, with neither specific reference to any particular 
method for constructing the solutions nor additional regularity assumptions. The Schwartz lemma 
infers from (39.1) that the distribution 

dtr^{u{t, x)) + Vj; • q(u(t, x)) 

is in fact a Radon measure, that is, the field {r]{u{t,x)),q{u{t,x))) is a divergence-measure field. 
Then there exists /x^ G A^(R+ x R") with /x^ < such that 

(39.2) diy^t,x)iviu{t,x)),q{u{t,x))) = /z,,. 

For any L°° entropy solution u, it is first indicated in Chen [Ch91] that, if the system is endowed 
with a strictly convex entropy, then, for any entropy pair {r],q), there exists /i^ G A1(R+ x R") 
such that 

(39.3) div(^t,x){v{u{t,x)),q{u{t,x))) = /i^. 
We introduce a functional on any oriented surface S: 

(39.4) MS)= I {Tl{u),q{u))-vd7e, 

Js 

where {rj{u), q{u))-v is the normal trace in the sense of Theorem 21.1, since {ri{u), q{u)) £ T>j\4]^^{M.+ X 
R"). It is easy to check that the functional defined by (39.4) is a Cauchy flux in the sense of 
Definition 32.2. 

39.2. Definition (Cauchy Entropy Fluxes). A functional defined by (39.4) is called a Cauchy 
entropy flux with respect to the entropy r/. 



40 GUI-QIANG CHEN MONICA TORRES WILLIAM P. ZIEMER 

In particular, when 77 is convex, then 

J'r^iS) > 

for any oriented surface S. Furthermore, we can reformulate the balance law of entropy from the 
recovery of an entropy production by capturing entropy dissipation. 

On the other hand, it is clear that understanding more properties of divergence-measure fields can 
advance our understanding of the behavior of entropy solutions for hyperbolic conservation laws and 
other related nonlinear equations by selecting appropriate entropy pairs. As examples, we refer the 
reader to [CF99a, CF99b, CF03, CW02] for the stability of Riemann solutions, which may contain 
rarefaction waves, contact discontinuities, and/or vacuum states, in the class of entropy solutions of 
the Euler equations for gas dynamics; to [CF99c, CW02] for the decay of periodic entropy solutions 
for hyperbolic conservation laws; to [CROO, VaOl] for the initial and boundary layer problems for hy- 
perbolic conservation laws; to [CF99a, CT05] for the initial-boundary value problems for hyperbolic 
conservation laws; and to [BFK02, MPT02] for nonlinear degenerate parabolic-hyperbolic equations. 

It is hoped that the theory of divergence-measure fields could be used to develop techniques in 
entropy methods, measure-theoretic analysis, partial differential equations, and related areas. 

Acknowledgments. The authors would like to thank Luis Caffarelli, Constantine Dafermos, Willi 

Jiiger, Fanghua Lin, Leon Simon, and David Swanson for stimulating and fruitful discussions. Gui- 
Qiang Chen's research was supported in part by the National Science Foundation under Grants 
DMS-0505473, DMS-0244473, and an Alexander von Humboldt Foundation Fellowship. Monica 
Torres's research was supported in part by the National Science Foundation under grant DMS- 
0540869. 



References 

L. Aiiibrosio, N. Fusco, and D. Pallara. Functions of Bounded Variation and Free Discontinuity Problems. 
Oxford Mathematical Monographs. The Clarendon Press, Oxford University Press: New York, 2000. 
G. Anzellotti. Pairings between measures and bounded functions and compensated compactness. Ann. 
Mat. Pura Appl. (4), 135:293-318 (1984), 1983. 

N. Aronszajn and K. T. Smith. Functional spaces and functional completion. Ann. Inst. Fourier. Grenoble, 
6:125-185, 1955-1956. 

C. Baiocchi and A. Capelo. Variational and Quasivariational Inequalities. John Wiley Sz Sons Inc.: New 
York, 1984. 

A. S. Besicovitch. On the definition of the area of a surface by means of inscribed polyhedra. J. London 
Math. Soc., 19:138-141, 1944. 

J. Bourgain and H. Brczis, On the equation divY = / and applications to control of phases, J. Amer. 

Math. Soc , 16(1): 393-426, 2002. 

F. Brezzi and M. Fortin. Mixed and Hybrid Finite Element Methods. Springer- Verlag: New York, 1991. 
Yu. D. Burago and V. G. Maa'ya. Potential Theory and Function Theory for Irregular Regions. Trans- 
lated from Russian. Seminars in Mathematics, V. A. Steklov Mathematical Institute, Leningrad, Vol. 3. 
Consultants Bureau, New York, 1969. 

R. Blirger, H. Prid, and K. H. Karlsen. On a free boundary problem for a strongly degenerate quasi- 
linear parabolic equation with an application to a model of pressure filtration. SIAM J. Math. Anal, 
34(3):611-635 (electronic), 2002. 

A. L. Cauchy. Rechcrches sur I'cquilibrc et Ic mouvement interieur des corps solides ou fiuides, elastiques 
or non elastiques. Bull. Soc. Philomathique, 10(2):9-13, 1823. 

A. L. Cauchy. Da la pression ou tension dans un corps solide. Exercises de Matehematiques, 2(2):42-56, 
1827. 

G. -Q. Chen. Hyperbolic systems of conservation laws with a symmetry. Commun. Partial Diff. Eqs., 
16:1461-1487, 1991. 

G.-Q. Chen and H. Prid. Divergence-measure fields and hyperbolic conservation laws. Arch. Ration. Mech. 
Anal, 147(2):89-118, 1999. 

G.-Q. Chen and H. Frid. Large-time behavior of entropy solutions of conservation laws. J. Diff. Eqs., 
152(2):308-357, 1999. 

G.-Q. Chen and H. Frid. Decay of entropy solutions of nonlinear conservation laws. Arch. Rational Mech. 
Anal, 146(2):95-127, 1999. 



[AFPOO] 

[Anz83] 

[AS56a] 

[BC84] 

[Bes44] 

[BB02] 

[BF91] 
[BM69] 

[BFK02] 

[Cal823] 
[Cal827] 
[Ch91] 
[CF99a] 

[CF99b] 
[CF99c] 



GAUSS-GREEN THEOREM, SETS OF FINITE PERIMETER, AND BALANCE LAWS 



41 



[CF03] G.-Q. Chen and H. Frid. Extended divergence-measure fields and the Euler equations for gas dynamics. 

Commun. Math. Phys., 236(2):251-280, 2003. 
[CROO] G.-Q. Ghen and M. Rascle. Initial layers and uniqueness of weaJc entropy solutions to hyperbolic conser- 
vation laws. Arch. Rational Mech. Anal, 153(3):205-220, 2000. 
[CT05] G.-Q. Chen and M. Torres. Divergence-measure fields, sets of finite perimeter, and conservation laws. Arch. 

Rational Mech. Anal, 175(2):245-267, 2005. 
[CTZb] G.-Q. Chen, M. Torres, and W. P. Zicrncr. Measure-theoretic analysis and nonlinear conservation laws. 

Pure Appl. Math. Quarterly, 2007 (to appear). 
[CW02] G.-Q. Chen and D.-H. Wang. The Cauchy problem for the Euler equations for compressible fluids. In 

Handbook of Mathematical Fluid Dynamics, Vol. I, pages 421-543. North-Holland: Amsterdam, 2002. 
[Daf98] C. M. Dafermos. Balance laws in continuum physics. In: Advanced Topics in Theoretical Fluid Mechanics 

(Paseky nad Jizerou, 1997), Pitman Res. Notes Math. Ser., 392, pages 89-117. Longman: Harlow, 1998. 
[Daf05] C. M. Dafermos. Hyperbolic Conservation Laws in Continuum Physics, 2nd Ed., Springer- Verlag: Berlin, 

2005. 

[DG61a] E. De Giorgi. Complementi alia teoria della misura (n — l)-dimensionale in uno spazio n-dimensionale. 

Scminario di Matematica della Scuola Normale Superiore di Pisa, 1960-61. Editrice Tecnico Scientifica, 
Pisa, 1961. 

[DG61b] E. Dc Giorgi. Frontiere orientate di misura minima. Scminario di Matematica della Scuola Normale 

Superiore di Pisa, 1960-61. Editrice Tecnico Scientifica, Pisa, 1961. 
[DMM99] M. Degiovanni, A. Marzocci, and A. Muscsti. Cauchy fiuxcs associated with tensor fields having divergence 

measure. Arch. Rational Mech. Anal, 147(3): 197-223, 1999. 
[DPP04a] T. De Pauw and W. F. Pfeffer. The Gauss-Green theorem and removable sets for PDEs in divergence 

form. Adv. Math., 183(1): 155-182, 2004. 
[EG92] L. C. Evans and R. F. Gariepy. Measure Theory and Fine Properties of Functions. Studies in Advanced 

Mathematics. CRC Press: Boca Raton, FL, 1992. 
[Fcd45] H. Fcdcrcr. The Gauss-Green theorem. Trans. Amer. Math. Soc, 58:44-76, 1945. 
[Fcd58] H. Fcdcrcr. A note on the Gauss-Green theorem. Proc. Amer. Math. Soc, 9:447-451, 1958. 
[Fed59] H. Fcdcrcr. Curvature measures. Trans. Amer. Math. Soc, 93:418-491, 1959. 
[Fed66] H. Fcdcrcr. Two theorems in geometric measure theory. Bull Amer. Math. Soc, 72:719, 1966. 
[Fed68b] H. Federer. Some properties of distributions whose partial derivatives axe representable by integration. 

Bull Amer. Math. Soc, 74:183-186, 1968. 
[Fed69] H. Federer. Geometric Measure Theory. Die Grundlehren der Mathematischen Wissenschaften, Band 153. 

Springer- Verlag New York Inc.: New York, 1969. 
[FF60] H. Federer and W. H. Fleming. Normal and integral currents. Ann. Math. (2), 72:458-520, 1960. 
[FZ73a] H. Federer and W. P. Ziemer. The Lebesgue set of a function whose distribution derivatives are p-th power 

summablc. Indiana Univ. Math. J., 22:139-158, 1972/73. 
[Fle60] W. H. Fleming. Functions whose partial derivatives are measures. Illinois J. Math., 4:452-478, 1960. 
[FL71] K. O. Friedrichs and P. D. Lax. Systems of conservation equations with a convex extension. Proc. Nat. 

Acad. Sci. U.S.A., 68:1686-1688, 1971. 
[Fug55] B. Puglede. On a theorem of F. Riesz. Math. Scand., 3:283-302 (1956), 1955. 

[Fug60] B. Puglede. Extremal Length and Closed Extensions of Partial Differential Operators. Jul. Gjellerups 

Boghandel, Copenhagen, 1960. 
[Giu77] E. Giusti. Minimal Surfaces and Functions of Bounded Variation. With notes by G. H. Williams, Notes 

on Pure Mathematics, 10, Department of Pure Mathematics, Australian National University, Canberra, 

1977. 

[GNN61] C. Goffman, C. J. Ncugebauer, and T. Nishiura. Density topology and approximate continuity. Duke Math. 
J., 28:497-505, 1961. 

[GW61] C. Goffman and D. Waterman. Approximately continuous transformations. Proc Amer. Math. Soc, 
12:116-121, 1961. 

[GM76] M. E. Gurtin and L. C. Martins. Cauchy's theorem in classical physics. Arch. Rational Mech. Anal, 
60(4):305-324, 1975/76. 

[GW67] M. E. Grurtin and W. O. Williams. An axiomatic foundation for continuum thermodynamics. Arch. Ra- 
tional Mech. Anal, 26:83-117, 1967. 

[GW71] M. E. Gurtin and W. O. Williams. On the first law of thermodynamics. Arch. Rational Mech. Anal, 
42:77-92, 1971. 

[GWZ86] M. E. Gurtin, W. O. Williams, and W. P. Zicrncr. Geometric measure theory and the axioms of continuum 

thermodynamics. Arch. Rational Mech. Anal, 92(l):l-22, 1986. 
[Gus60] W. Gustin. Boxing inequalities. J. Math. Mech., 9:229-239, 1960. 
[He07] J. Heinonen. Nonsmooth calculus. Bull Amer. Math. Soc, 1:632-232, 2007. 

[JN90] W. B. Jurkat and D. J. F. Nonnenmacher. The general form of Green's theorem. Proc. Amer. Math. Soc, 
109(4): 1003-1009, 1990. 



42 



GUI-QIANG CHEN 



MONICA TORRES 



WILLIAM P. ZIEMER 



[JN94a] W. B. Jurkat and D. J. F. Nonnenmacher. An axiomatic theory of non-absolutely convergent integrals in 
R". Fund. Math., 145(3):221-242, 1994. 

[JN94b] W. B. Jurkat and D. J. F. Nonnenmacher. A generahzed n-dimensional Riemann integral and the diver- 
gence theorem with singularities. Acta Sci. Math. (Szeged), 59(l-2):241-256, 1994. 

[JN94c] W. B. Jurkat and D. J. F. Nonnenmacher. A theory of non-absolutely convergent integrals in R" with 
singularities on a regular boundary. Fund. Math., 146(1) :69-84, 1994. 

[JN95] W. B. Jurkat and D. J. F. Noniicmriaclicr. A Hake-type property for the t'l-integral and its relation to 
other integration processes. Czechoslovak Math. J., 45(120){3):465-472, 1995. 

[JNm95] W. B. Jurkat and D. J. F. Nonnenmacher. The fundamental theorem for the j^i-integral on more general 
sets and a corresponding divergence theorem with singularities. Czech. Math. J., 45(120)(l):69-77, 1995. 

[Lew88a] J. L. Lewis. Uniformly fat sets. Trans. Amer. Math. Sac, 308(1):177-196, 1988. 

[LY02] F. H. Lin and X. P. Yang. Geometric Measure Theory — An Introduction. Science Press; Beijing; Interna- 
tional Press: Boston, MA, 2002. 

[LRX07] A. Lim, B. Rodrigues, and Zhou Xu. A m-parallel crane scheduling problem with a non-crossing constraint. 
Naval Res. Logist., 54(2):115-127, 2007. 

[MPT02] C. Mascia, A. Porretta, and A. Terracina. Nonhomogencous Dirichlct problems for degenerate parabolic- 
hyperbolic equations. Arch. Rational Mech. Anal, 163(2):87-124, 2002. 

[Mor47] A. P. Morse. Perfect blankets. Trans. Amer. Math. Soc, 61:418-442, 1947. 

[Nol59] W. Noll. The foundations of classical mechanics in the light of recent advances in continuum mechanics. In: 
The Axiomatic Method. With Special Reference to Geometry and Physics. Proceedings of an International 
Symposium held at the Univ. of Calif., Berkeley, Dec. 26, 1957- Jan. 4, 1958 (edited by L. Henkin, P. 
Suppes and A. Tarski), Studies in Logic and the Foundations of Mathematics, pages 266-281, North- 
Holland Publishing Co.: Amsterdam, 1959. 

[Nol72] W. Noll. A new mathematical theory of simple materials. Arch. Rational Mech. Anal., 48:1-50, 1972. 

[Nol73] W. Noll. Lectures on the foundations of continuum mechanics and thermodynamics. Arch. Rational Mech. 
Anal, 52:62-92, 1973. 

[N0I86] W. Noll. Continuum mechanics and geometric integration theory. In: Categories in Continuum Physics 
(Buffalo, N.Y., 1982), Lecture Notes in Math., 1174, pages 17-29. Springer: Berlin, 1986. 

[Nol87] W. Noll. Finite- Dimensional Spaces. Vol. I, Mechanics: Analysis, 10. Maxtinus Nijhoff Publishers: Dor- 
drecht, 1987. 

[Non95] D. J. F. Nonnenmacher. Sets of finite perimeter and the Gauss-Green theorem with singularities. J. London 

Math. Soc. (2), 52(2):335-344, 1995. 
[Pfe] W. F. Pfeffer. Distributions for which divv = f has a continuous solution. To Appear. 

[PfeOl] W. F. Pfeffer. The Stokes theorem for the generalized Riemann integral. Real Anal. Exchange, 26(2):623- 

636, 2000/01. 

[Pfe90] W. F. Pfeffer. Divergence theorem for vector fields with singularities. In: New Integrals (Coleraine, 1988), 

Lecture Notes in Math., 1419, pp. 150-166, Springer: Berlin, 1990. 
[Pfe05a] W. F. Pfeffer. The Gauss-Green theorem in the context of Lebesgue integration. Bull. London Math. Soc, 

37(l):81-94, 2005. 

[Sch07] F. Schuricht. A new mathematical foundation for contact interactions in continuum physics. Arch. Rational 

Mech. Anal, 184: 495-551, 2007. 
[Sim83] L. Simon. Lectures on Geometric Measure Theory, Proceedings of the Centre for Mathematical Analysis, 

3, Australian National University, Canberra, 1983. 

[Ste93] E. M. Stein. Harmonic Analysis: Real-Variable Methods, Orthogonality, and Oscillatory Integrals. Prince- 
ton Mathematical Series, 43. Princeton University Press: Princeton, NJ, 1993. 

[VaOl] A. Vasseur. Strong traces for solutions of multidimensional scalar conservation laws. Arch. Rational Mech. 
Anal, 160(3):181-193, 2001. 

[Vol67a] A. I. Vol'pert. Spaces BV and quasilinear equations. Mat. Sb. (N.S.), 73 (115):255-302, 1967. 

[VH69a] A. I. Vol'pert and S. I. Hudjaev. The Cauchy problem for second order quasilinear degenerate parabolic 
equations. Mat. Sb. (N.S.), 78 (120):374-396, 1969. 

[Whi57] H. Whitney. Geometric Integration Theory. Princeton University Press: Princeton, N. J., 1957. 

[Zic67] W. P. Ziemer. Extremal length and conformal capacity. Trans. Amer. Math. Soc, 126:460—473, 1967. 

[Zic69] W. P. Ziemer. Extremal length and p-capacity. Michigan Math. J., 16:43—51, 1969. 

[Zie70] W. P. Ziemer. Extremal length as a capacity. Michigan Math. J., 17:117-128, 1970. 

[Zie83] W. P. Ziemer. Cauchy flux and sets of finite perimeter. Arch. Rational Mech. Anal, 84(3):189-201, 1983. 
[Zie89] W. P. Ziemer. Weakly Differentiable Functions, Graduate Texts in Mathematics, 120, Springer- Verlag: 
New York, 1989. 



(C.-Q. Chen) Department of Mathematics, Northwestern University, 2033 Sheridan Road, Evanston, IL 
60208-2730, USA. http://www.math. northwestern. edu/"gqchen 
E-mail address: gqchenamath.northwestern.edu 



GAUSS-GREEN THEOREM, SETS OF FINITE PERIMETER, AND BALANCE LAWS 



43 



(M. Torres) Department of Mathematics, Purdue University, 150 N. University Street, West Layayette, 
IN 47907-2067, USA. http://www.math.purdue.edu/-torres 
E-mail address: torres@math.purdue.edu 

(W. Ziemer) Department of Mathematics, Indiana University, Rawles Hall, Bloomington, IN 47405, 
USA. HTTP : / / www .indiand.edu/"'ziemer 
E-mail address: ziemerQindiana.edu 



